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Abstract 

We perform a dedicated study of the four-fermion production process e~e + — > 
fi^u^udX near the W pair-production threshold in view of the importance of 
this process for a precise measurement of the W boson mass. Accurate theoret- 
ical predictions for this process require a systematic treatment of finite-width 
effects. We use unstable-particle effective field theory (EFT) to perform an 
expansion in the coupling constants, Tw/Mw, and the non-relativistic velocity 
v of the W boson up to next-to-leading order in Tw/Mw ~ ct ew ~ v 2 . We 
find that the dominant theoretical uncertainty in My/ is currently due to an 
incomplete treatment of initial-state radiation. The remaining uncertainty of 
the NLO EFT calculation translates into 5My/ « 10 - 15 MeV, and to about 
5 MeV with additional input from the NLO four-fermion calculation in the full 
theory. 



1 Introduction 



The mass of the W gauge boson is a key observable in the search for virtual-particle 
effects through electroweak precision measurements. Its current value, Mw = (80.403 ± 
0.029) GeV [1], is determined from a combination of continuum W pair-production at 
LEPII and single- W production at the Tevatron. 1 Further measurements of single- W 
production at the LHC should reduce the error by a factor of two. Beyond LHC it has 
been estimated that an error of 6 MeV could be achieved by operating an e"e + collider in 
the vicinity of the W pair-production threshold [3]. This estimate is based on statistics 
and the performance of a future linear collider, and it assumes that the cross section is 
known theoretically to sufficient accuracy so that its measurement can be converted into 
one of Mw- In reality, achieving this accuracy is a difficult theoretical task, requiring 
the calculation of loop and radiative corrections. Since the W bosons decay rapidly, 
this calculation should be done for a final state of sufficiently long-lived particles, rather 
than for on-shell W pair-production. A systematic treatment of finite-width effects is 
therefore needed. 

In this paper we investigate in detail the inclusive four-fermion production process 

e~(pi) e + (j» 2 ) -> H~ v^ud + X (1) 

in the vicinity of the W pair-production threshold, i.e. for s = {P1+P2) 2 ~ ^M w . Here X 
denotes an arbitrary flavour-singlet state (nothing, photons, gluons, ...). No kinematic 
cuts shall be applied to the final state. In this kinematical regime the process (1) is 
primarily mediated by the production of two resonant, non-relativistic W bosons with 
virtuality of order 

k 2 — M w ~ M w v 2 ~ M W T W < M w , (2) 

one of which decays into leptons, the other into hadrons. Here we have introduced 
the non-relativistic velocity v, and the W decay width Yw- We perform a systematic 
expansion of the total cross section in the small quantities 



S-4:M W 2 L W 



corresponding to a (re-organized) loop expansion and a kinematic expansion. All three 
expansion parameters are of the same order, and for power-counting purposes we denote 
them collectively as 5. Our calculation is accurate at next-to- leading order (NLO). 
Note that resonant processes such as (1) are complicated by the need to account for 
the width of the intermediate unstable particles to avoid kinematic singularities in their 
propagators. The expansion in the electroweak coupling a ew = a.js 2 w is therefore not 
a standard loop expansion, (a denotes the electromagnetic coupling, and = sin 2 6 W 
with 6 W the Weinberg angle.) 



1 This value refers to the definition of the W mass from a Breit-Wigner parameterization with a 
running width as it is adopted in the experimental analyses. It is related to the pole mass Mw used in 
this paper by [2] M w - M w = T W /(2M W ) + 0(a 3 ew ). 
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NLO calculations of four-fermion production have been done already some time ago 
in the continuum (not near threshold) in the double-pole approximation for the two W 
propagators [4-6] or with further simplifications [7, 8] . This approximation was supposed 
to break down for kinematic reasons in the threshold region. Thus, when this project 
was begun [9], there existed only LO calculations in the threshold region as well as 
studies of the effect of Coulomb photon exchanges [10,11], rendering the effective field 
theory approach [12-14] the method of choice for the NLO calculation. Meanwhile a full 
NLO calculation of four-fermion production has been performed in the complex mass 
scheme [15,16] without any kinematic approximations, and for the fully differential cross 
sections in the continuum or near threshold. This is a difficult calculation that required 
new methods for the numerical evaluation of one-loop six-point tensor integrals. In 
comparison, our approach is computationally simple, resulting in an almost analytic 
representation of the result. The drawback is that our approach is not easily extended 
to differential cross sections. Nevertheless, we believe that a completely independent 
calculation of NLO four-fermion production is useful, and we shall compare our result 
to [15] in some detail. Having a compact analytic result at hand is also useful for an 
investigation of theoretical uncertainties. Note that while the full four-fermion NLO 
calculation [15, 16] is a priori of the same accuracy in Tw/Mw as the NLO effective- 
theory result, it includes a subset of higher order terms in the EFT expansion. We 
discuss the relevance of these higher order terms at the end of this paper. 

The organization of the paper is as follows. In Section 2 we explain our method of 
calculation. We focus on aspects of unstable-particle effective theory that are specific to 
pair production near threshold and refer to [13] for those, which are in complete analogy 
with the line-shape calculation of a single resonance. The section ends with a list of all 
terms that contribute to the NLO result. We construct the effective-theory expansion of 
the tree approximation to the four-fermion cross section in Section 3. Of course, this cal- 
culation can be done nearly automatically without any expansions with programs such 
as Whizard [17], CompHep [18,19] or MadEvent [20,21]. The purpose of this section is 
to demonstrate the convergence of the expansion towards the "exact" tree-level result, 
and to provide analytic expressions for those terms that form part of the NLO calcula- 
tion near threshold. In Section 4 we calculate the radiative corrections required at NLO. 
These consist of hard loop corrections to W pair-production and W decay, of Coulomb 
corrections up to two photon exchanges, and soft-photon corrections. The entire calcu- 
lation is done setting the light fermion masses to zero, which is a good approximation 
except for the initial-state electrons, whose mass is relevant, since the cross section is not 
infrared-safe otherwise. In Section 5 we describe how to transform from the massless, 
"partonic" cross section to the physical cross section with finite electron mass, includ- 
ing a resummation of large logarithms \n(s/m%) from initial-state radiation. Assembling 
the different pieces we obtain the full inclusive NLO four-fermion cross section in terms 
of compact analytic and numerical expressions. In Section 6 we perform a numerical 
evaluation of the NLO cross section, estimate the final accuracy, and compare our result 
to [15], obtaining very good agreement. We find that the dominant theoretical uncer- 
tainty in Mw is currently due to an incomplete treatment of initial-state radiation. The 
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remaining uncertainty of the NLO EFT calculation translates into SM W pz 10 - 15 MeV, 
and to about 5 MeV with additional input from the NLO four-fermion calculation in the 
full theory. We conclude in Section 7. Some of the lengthier equations are separated 
from the main text and provided in Appendices A and B. 



We extract the inclusive cross section of the process (1) from the appropriate cuts of 
the e~e + forward-scattering amplitude. For inclusive observables, where one integrates 
over the virtualities of the intermediate resonances, the propagator singularity poses 
no difficulty, if the integration contours can be deformed sufficiently far away from the 
singularity. This is not possible, however, for the calculation of the line-shape of a single 
resonance, and for pair production near threshold (the pair production equivalent of the 
resonance region), where the kinematics does not allow this deformation. The width of 
the resonance becomes a relevant scale, and it may be useful to separate the dynamics 
at this scale from the dynamics of the short distance fluctuations at the scale of the 
resonance mass by constructing an effective field theory. 

2.1 Unstable-particle effective theory for pair production near 
threshold 

The following formalism resembles rather closely the formalism described in [12,13]. The 
generalization from a scalar to a vector boson resonance is straightforward. The pair- 
production threshold kinematics implies a change in power counting that is analogous 
to the difference between heavy-quark effective theory and non-relativistic QCD. 

In W pair-production the short-distance fluctuations are given by hard modes, whose 
momentum components are all of order Myy. After integrating out the hard modes, the 
forward-scattering amplitude is given by [13] 



iA = J2 ^(e~e + |T[^t(o),o(O( x )]| e - e +) + ^( e - e +|, W(o)| e -e + ). (4) 



The operators O p '(x) (O p (x)) in the first term on the right-hand side produce (destroy) 
a pair of non-relativistic W bosons. The second term accounts for the remaining non- 
resonant contributions. The matrix elements are to be computed with the effective 
Lagrangian discussed below and the operators include short-distance coefficients due to 
the hard fluctuations. Note that there is no separate term for production of one resonant 
and one off-shell W, since for such configurations the integrations are not trapped near 
the singularity of the W propagator. These configurations are effectively short-distance 
and included in the non- resonant production-decay operators O^J(0). 

The effective Lagrangian describes the propagation and interactions of two non- 
relativistic, spin-1 fields Q l ± representing the nearly on-shell (potential) W ± modes; two 
sets of collinear fields for the incoming electron and positron, respectively; and potential 
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and collinear photon fields. The corresponding momentum scalings in the center-of-mass 
frame are: 

potential (p) : k ~ M w 5, \k\ ~ M W V6 

soft (s) : k ~ \k\ ~ M w <5 (5) 

collinear (c) : &o ~ , /c 2 ~ M^S. 

The small parameter <5 is either the non-relativistic velocity squared, v 2 , related to 
(s — 4M^)/(4M^), or Tw/Mw ~ « e », since the characteristic virtuality is never para- 
metrically smaller than M W T W for an unstable W. The interactions of the collinear 
modes are given by soft-collinear effective theory [22-24]. There is nothing specifically 
new related to collinear modes in pair production, and we refer to [13] for further de- 
tails. As far as the next-to-leading order calculation is concerned, the soft-collinear 
Lagrangian allows us to perform the standard eikonal approximation for the interaction 
of soft photons with the energetic electron (positron) in the soft one-loop correction. 

The Lagrangian for the resonance fields is given by the non-relativistic Lagrangian, 
generalized to account for the instability [9,25]. The terms relevant at NLO are 



'nrqed — ^ 



" 2M W 2 ) a a 8M^ 



(6) 



Here Q l + and Q l _ (i = 1,2,3) are non-relativistic, spin-1 destruction fields for particles 
with electric charge ±1, respectively. The interactions with photons is incorporated 
through the covariant derivative D^Vt % ± = (d^ =F ieA^)VL l ± . The effective theory does not 
contain fields for the other heavy particles in the Standard Model, the Z and Higgs 
bosons, and the top quark. Their propagators are always off-shell by amounts of order 
M^y and therefore their effect is encoded in the short-distance matching coefficients. In 
a general i?g-gauge this also applies to the pseudo-Goldstone (unphysical Higgs) fields, 
except in 't Hooft-Feynman gauge £ = 1, where the scalar W and unphysical charged 
pseudo-Goldstone modes have masses My/ and can also be resonant. However, the two 
degrees of freedom cancel each other, leading to the same Lagrangian (6) describing the 
three polarization states of a massive spin-1 particle. The effective Lagrangian has only 
a U(l) electromagnetic gauge symmetry as should be expected at scales far below M\y. 
However, since the short-distance coefficients of the Lagrangian and all other operators 
are determined by fixed-order matching of on-shell matrix elements to the full Standard 
Model, they are independent of the gauge parameter in i?^-gauge by construction. The 
often quoted gauge-invariance problems in the treatment of unstable particles arise only 
if one performs resummations of perturbation theory in gauge-dependent quantities such 
as propagators. 

The matching coefficient A in (6) is obtained from the on-shell two-point function of 
a transverse W boson. "On-shell" here refers to the complex pole determined from 



s-M 2 w -I^(s)=Q (7) 
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with M\y any renormalized mass parameter, and U^(q 2 ) the renormalized, transverse 
self-energy. The solution to this equation, 

s = M w - %M W Y W) (8) 

defines the pole mass and the pole width of the W. The matching coefficient is then 
given by 

S — M w pole scheme 

= — - = —li-w- i y ) 

M w 

In the remainder of the paper, we adopt a renormalization convention where M w is the 
pole mass M w , in which case A is purely imaginary. With D° ~ M w 5, D 2 ~ M w 5, and 
A ~ My/5, we see that the first bilinear term in (6) consists of leading-order operators, 
while the second is suppressed by one factor of 5, and can be regarded as a perturbation. 
Accordingly, the propagator of the fl± fields is 

(10) 



k° 



k 2 _ A 
2M W 2 



The effective theory naturally leads to a fixed-width form of the resonance propaga- 
tor. Note that it would be sufficient to keep only the one-loop expression for A in the 
propagator, and to include higher-order corrections perturbatively. 

Loop diagrams calculated using the Lagrangian (6) receive contributions from soft 
and potential photons. 2 Since the potential photons do not correspond to on-shell parti- 
cles, they can be integrated out, resulting in a non-local (Coulomb) potential, analogous 
to potential non-relativistic QED [27]. Up to NLO the required PNRQED Lagrangian 
is 



-PNRQED — 



fit* ( iD° + — — W + ( * ~ MWA)2 



(11) 



+ J d 3 f [ife: ] ( x + f ) (- r 

Only the (multipole-expanded) soft photon A° s (t, 0) appears in the covariant derivative 
Dg. The potential W field has support in a region ~ 5 _1 in the time direction and in 
a region ~ 5~ x / 2 in each space direction, hence the measure d 4 x in the action scales as 
5 _5//2 . Together with d ~ 5 we find from the kinetic term that il l T ~ 5 3 / 4 . Analogously 
we find that the non-local Coulomb potential scales as a/\^5 ~ a/v. Since we count 
a ~ v 2 , the Coulomb potential is suppressed by v, or a 1 / 2 , and need not be resummed, 
in contrast to the case of top-quark pair-production near threshold. However, with this 
counting the Coulomb enhancement introduces an expansion in half-integer powers of 
the electromagnetic coupling, the one-loop Coulomb correction being a "N 1 / 2 LO" term. 



2 What we call "soft" here, is usually termed "ultrasoft" in the literature on non-relativistic QCD. 
There are further modes (called "soft" there) with momentum k ~ M\y\/S [26]. In the present context 
these modes cause, for instance, a small modification of the QED Coulomb potential due to the one-loop 
photon self-energy, but these effects are beyond NLO. 
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Figure 1: Diagrams contributing to the tree- level matching of O 



(0) 



2.2 Production vertex, production-decay vertices and the lead- 
ing-order cross section 

We now turn to the production and production-decay operators appearing in the rep- 
resentation (4) of the forward-scattering amplitude. The lowest-dimension production 
operator must have field content (e C2 e cl ) (fPfi+), where the subscripts on the electron 
fields stand for the two different direction labels of the collinear fields. The short-distance 
coefficients follow from matching the expansion of the renormalized on-shell matrix el- 
ements for e~e + — > W~W + in the small relative W momentum to the desired order 
in ordinary weak-coupling perturbation theory. The on-shell condition for the W lines 



implies that their momentum satisfies k[ 



M^r + M\yA, but in a perturbative 



matching calculation this condition must be fulfilled only to the appropriate order in a 
and S. On the effective-theory side of the matching equation one also has to add a factor 
^/2M^w~ 1/2 with 

1 + M «t +P T (12) 



zu 



M 2 W 



for each external Q line [13]. 3 At tree-level, and at leading order in 5, zu^ 1 = 1. 

Thus we are led to consider the tree-level, on-shell W pair-production amplitude 
shown in Figure 1. To leading order in the non-relativistic expansion the s-channel 
diagrams vanish and only the helicity configuration e^e\ contributes. The corresponding 
operator (including its tree-level coefficient function) reads 



O (o) 



M 2 W 



(13) 



where we have introduced the notation = a l b> + a 3 b l and the unit-vector n for the 
direction of the incoming electron three- momentum pf. For completeness we note that 
the emission of collinear photons from the W or collinear fields of some other direction, 
which leads to off-shell propagators, can be incorporated by adding Wilson lines to the 
collinear fields, implying the form (e C2j i,W C2 7' l n^W ( | 1 e Cl) £,). However, these Wilson lines 
will not be needed for our NLO calculation, since the collinear loop integrals vanish (see, 
however, Section 5). 



3 This is the well-known (E/M) 1 ^ 2 factor, which accounts for the normalization of non-relativistic 
fields, generalized to unstable particles and general mass renormalization conventions. 
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Figure 2: Leading-order effective-theory diagram for the forward-scattering amplitude. 



The leading contribution from the potential region to the forward-scattering ampli- 
tude is given by the expression 

MS = / d 4 x (ele + R \T[ i O^\0) i O^\x)]\e- L e + R ). (14) 

This corresponds to the one- loop diagram shown in Figure 2, computed with the ver- 
tex (13) and the propagator (10). We can use power counting to estimate the magnitude 
of the leading-order amplitude prior to its calculation. With e CijL ~ 5 1/2 , ~ <5 3/4 
the production operator scales as ~ a5 5 ^ 2 . The integration measure scales as 
J d A x ~ <5~ 5 / 2 in the potential region and the external collinear states are normalized 

as |e T ) ~ hence A^ R ~ a 2 5 1 ^ 2 . This expectation is confirmed by the explicit 

calculation of the one-loop diagram: 

= ^ (P2 - I^W-XPi - Wh 3] \P2-) 

r d d r 1 
X J {2n) d ( o _ Z a \( E _ o _ ZEL - 

y 2M W 2 J \y 1 2M W 2 J 

= ~ Ai7iaL (15) 

Here we have defined E = y/s — 2Mw- We adopted the standard helicity notation 
\p±) = ^r~u(p), and used A = — iT w , valid in the pole scheme, in the last line. 
The fermion energies are set to M w in the external spinors. The calculation has been 
performed by first evaluating the r° integral using Cauchy's theorem, and the trace 
(P2 — \n)- l ^\pi— )(pi — |n' l 7^|p2 — ) = 16(1 — e)My V . The remaining \ f\ integral contains 
a linear divergence that is, however, rendered finite by dimensional regularization (with 
d — 4 — 2e) so the d — > 4 limit can be taken. The numerical comparison of (15) to the 
full tree-level result and the convergence of the effective-theory approximation will be 
discussed in Section 3. 

Taking the imaginary part of (15) does not yield the cross section of the four-fermion 
production process (1) with its flavour-specific final state. At leading order the cor- 
rect result is given by multiplying the imaginary part with the leading-order branching 
fraction product Bt^°\W~ — > fi^u^)Br^(W + — > ud) = 1/27. This procedure can be 
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Figure 3: Cut one- loop diagrams contributing to non- resonant production-decay operator 
matching. 



justified as follows. The imaginary part of the non-relativistic propagator obtained by 
cutting an Q line is given by 

l r {0) 12 



2M W ^2 ^ 2M W / ^ 



The propagator of the Q± line implicitly includes a string of self-energy insertions. Taking 
the imaginary part amounts to performing all possible cuts of the self-energy insertions 
while the unstable particle is not cut [28]. To obtain the total cross section for a flavour- 
specific four-fermion final state, only the cuts through these specific fermion lines have 
to be taken into account. At the leading order this amounts to replacing T$ in the 
numerator of (16) by the corresponding partial width, here and T^, respectively, 

while the total width is retained in the denominator. The leading-order cross section is 
therefore 



M _ 1 . „(0) 4tt« 2 



LR 27s LR 274« 



M w 

(0) 



(17) 



The unpolarized cross section is given by (x^/4, since the other three helicity combina- 
tions vanish. 

The leading contribution from non-resonant pro duct ion- decay operators O^J to (4) 
arises from four-electron operators of the form 



° ( £ = T?H e ~ C i r i e C2 )(^r 2 e cl ), (18) 



Ml 



where r 1; T 2 are Dirac matrices. If Cj^ ~ a n , the contribution to the forward-scattering 
amplitude scales as at n . This should be compared to A^ R ~ a 2 5 1 ^ 2 . The calculation 
of the short-distance coefficients Cj^ is performed in standard fixed-order perturbation 
theory in the full electroweak theory. The W propagator is the free propagator, since the 
self-energy insertions are treated perturbatively. The leading contribution to the forward- 
scattering amplitude arises from the one-loop diagrams shown in Figure 3. We will 
calculate the imaginary part of the short- distance coefficients cj£? by evaluating the cut 
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diagrams. The calculation of cuts corresponding to tree amplitudes is most conveniently 
performed in unitary gauge with W propagator — i(g IJ , l/ — k^k^/M^) / (k 2 — + ie). To 
leading order in the expansion in 5, the cut one-loop diagrams in Figure 3 correspond 
to the production cross section of two on-shell W bosons directly at threshold, which 
vanishes. In fact, from an explicit representation of these one-loop diagrams it can be 
seen that the imaginary parts from the hard region vanish in dimensional regularization 
to all orders in the 5 expansion. Thus the leading imaginary parts of Cf} arise from 
two-loop diagrams of order a 3 . Just as the Coulomb correction the leading non- resonant 
(hard) contribution provides another N 1//2 LO correction relative to (15). 

2.3 Classification of corrections up to NLO 

We now give an overview of the contributions to the four-fermion cross section at N^LO 
and NLO. These consist of the short-distance coefficients of the non-relativistic La- 
grangian (11), of the production operators Op k \ and the four-electron operators Of} on 
the one hand; and corrections that arise in calculating the matrix elements in (4) within 
the effective theory on the other. 

2.3.1 Short-distance coefficients in the effective Lagrangian 

The effective Lagrangian (11) is already complete to NLO. The only non-trivial matching 
coefficient is A, which follows from the location of the W pole, which in turn can be 
computed from the expansion of the self-energy [13]. In the pole scheme, we require 
the NLO correction to the decay width Tw, defined as the imaginary part of the pole 
location, see (8), (9). At leading order, A^ = — iF$ with 4 

T% ] = - A a ew M w . (19) 
There are electroweak as well as QCD corrections to the W self-energy. We shall 

1 /2 

count the strong coupling a s as a ew ■ Thus the mixed Q CD-electro weak two-loop self- 
energy provides a N 1 / 2 LO correction to A, while at NLO we need the self-energy at 
orders o? ew and a ew a 2 . The QCD effects are included by multiplying the leading-order 
hadronic partial decay widths by the universal QCD correction for massless quarks [29], 

<Jqcd = 1 + — + 1.409-|, (20) 

with a s = a s (M w ) in the MS scheme. The electroweak correction to the pole-scheme 
decay width is denoted by T^ ew \ The explicit expression is given in Section 4.1. We 
therefore have 

A (3/ 2 ) = _ ir (V2) = _ t 2^ r (o ); A(2) = _ ?r « = _ % 

Oil 

(21) 

4 Here the masses of the light fermions are neglected, and the CKM matrix has been set to the unit 
matrix. 



r& ew) + 1.409 ^r (0) 
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These results refer to the total width, which appears in the propagator and the forward- 
scattering amplitude. The extraction of the flavour-specific process e~e + — > fi~u^udX 
will be discussed in Section 3.2. 



2.3.2 Matching coefficients of the production operators 

There are two sorts of corrections related to production operators: higher-dimensional 
operators suppressed by powers of 5, and one-loop corrections to the operators of lowest 
dimension such as (13). 

The higher-dimension production operators are of the form 

Of = -^W^L/RTHn, D)e L/R )(tTlg(D)n% (22) 

where V is some combination of Dirac matrices and JF and Q are functions of the covariant 
derivative D acting on the fields. (Here and below, we drop the collinear direction label 
on the electron fields, whenever they are obvious.) The short-distance coefficients of 
these operators follow from the expansion of appropriate on-shell amplitudes around the 
threshold. The expansion parameter is v ~ 5 1 / 2 . However, for the inclusive cross section 
there is no interference of the f-suppressed operator with the leading one, hence the 
correction from higher- dimension operators begins at NLO. Full results for the tree-level 
matching of the N 1,/2 LO production operators are given in [9]. The NLO contribution to 
the inclusive cross section is computed in Section 3.1. 

The one- loop correction to the matching coefficient of the production vertex (13) and 
the related operator with right-handed electrons requires to calculate the renormalized 
scattering amplitudes for e^e^ — > W + W~ and e R e^ — > W + W~ to NLO in ordinary 
weak coupling perturbation theory for the momentum configuration (pi + p 2 ) 2 = 4M^, 
i.e. directly at threshold. This generates the NLO production operator 

Of = ^jr [Cfn (e L l^e L ) + C^ L (e Rl We R )] . (23) 

The calculation of the coefficients Cf LR , Cf RL is discussed in Section 4.1. Note, however, 
that the one-loop correction Cf RL does in fact not contribute to the NLO cross section, 
since there is no leading-order contribution from the e R e^ helicity initial state, and no 
interference between LR and RL configurations. 



2.3.3 Matching coefficients of four-electron operators 

As discussed above the leading contributions from the non-resonant production-decay 
operators to the imaginary part of the forward scattering amplitude arise at N 1 / 2 LO, 
where the half-integer scaling arises from the absence of the threshold suppression v ~ 
^1/2 p resen t j n the LO cross section. The calculation of the cut 2-loop diagrams amounts 
to the calculation of the squared and phase-space integrated matrix element of the on- 
shell processes e~e + — > W~ud and e~e + — * jj,~u^W + in ordinary perturbation theory 
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(no "resummations" in internal W propagators). This includes contributions of what 
is usually called double-resonant (or CC03) diagrams, where one of the W propagators 
is in fact off-shell, as well as genuine single-resonant processes. In the terminology of 
the method of regions, these corrections are given by the hard-hard part of the two-loop 
forward-scattering amplitude. Since they contain all diagrams contributing to the tree- 
level scattering processes e"e + — > fi~u fl W + and e"e + — > W~ud, the matching coefficients 
are gauge invariant. Since only one W line is cut in the N 1//2 LO contributions, they can 
be viewed as systematic corrections to the narrow- width approximation. This calculation 
is performed in Section 3.3. 

To NLO in the power counting a 2 s ~ a ew we would have to compute also the NLO 
QCD corrections to e~e + — > W~ud (+<?). The corrections to the "double-resonant" 
(CC03) diagrams can be taken into account approximately by multiplying them with 
the one-loop QCD correction to the hadronic decay width. The corrections to the single- 
resonant diagrams require the full calculation. However, we shall find that the contribu- 
tion of the single- resonant diagrams to e~e + — > W~ud is numerically already small, so 
we neglect the QCD corrections. 

2.3.4 Calculations in the effective theory 

One-loop diagrams with insertions of subleading operators. The contributions in this 
class arise from evaluating the first term in (4) at one loop, see Figure 2, but with one 
insertion of the subleading bilinear terms in the Lagrangian (11), which correspond to 
kinetic energy and width corrections, or with production operator products Op^OjP 
and O^ 2 ^ 0^ 2 \ where is either a higher- dimension operator (22) or the one-loop 
correction (23). As already mentioned the N 1//2 LO products O^O^ 2 ^ vanish after 
performing the angular integrals. In the calculation discussed further in Section 3 we 
actually follow a different approach and directly expand the spin-averaged squared matrix 
elements rather than the amplitude before squaring, which would yield the individual 
production vertices. 

Coulomb corrections. A single insertion of the Coulomb potential interaction in the 
Lagrangian (11) contributes at N^LO. To NLO one has to calculate the double insertion 
into the leading-order amplitude from O^of^ and a single insertion into 0^0^ 2 \ The 
latter vanishes for the total cross section. There is no coupling of the potential photons 
to the collinear electrons and positrons, so there are no Coulomb corrections to the 
four-fermion operators. The Coulomb corrections are given in Section 4.2. 

NLO corrections from soft and collinear photons. To NLO one has to calculate two- 
loop diagrams in the effective theory arising from the coupling of the collinear modes 
and the potential W bosons to the soft and collinear photons contained in the NRQED 
Lagrangian (6) and the SCET Lagrangian. The cuts correspond to one-loop virtual 
and bremsstrahlung corrections to the leading-order cross section. In the terminology of 
the method of regions these are contributions from the soft-potential, the Ci-potential 
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and the C2-potential regions. They correspond to "non-factorizable corrections" and are 
discussed in Section 4.3. 



3 Expansion of the Born cross section 

This section serves two purposes. First, we calculate all NLO corrections to four-fermion 
production in the effective theory (EFT) except those related to loop corrections, which 
will be added in Section 4. Second, we investigate the convergence of the successive EFT 
approximations to what is usually referred to as the Born four-fermion production cross 
section. The two calculations are not exactly the same, since the implementation of the 
W width in the Born cross section is not unique. We define the "exact" Born cross section 
by the ten tree diagrams for e~e + — > fi~i?^ud, where the W propagators are supplied 
with a fixed-width prescription. The EFT calculation is done by expanding directly 
the forward-scattering amplitude. The relevant loop momentum regions are either all 
hard, or hard and potential. In the latter regions the two W propagators and the W 
interactions are described by the non-relativistic Lagrangian. The all-hard contributions 
correspond to the matching and matrix element of the four-electron operators. 



3.1 Expansion in the potential region 

We first reconsider the one-loop diagrams (before cutting) shown in Figure 3, where the 
loop momentum is now assumed to be in the potential region. The forward-scattering 
amplitude corresponding to these diagrams may be written as 

lA = I -0y!^ r ) P (h)P(k 2 ) : (24) 

where E = y/s-2M w , k x = M w v+r, k 2 = P-M w v-r, with v* 1 = (1,0) and P = pi+p 2 
the sum of the initial-state momenta. Here <&(E, r) is the square of the off-shell W pair- 
production amplitude at tree level, including the numerator {—g^ + k^k u /k 2 ) from the 
W propagators, and 

m = w^m^wW) (25) 

is the full renormalized (transverse) W propagator. 5 Writing the amplitude in the full 
theory with a resummed propagator is contrary to the spirit of effective field theory 
calculations, where the matching coefficients are obtained by fixed-order calculations. 
However, this allows us to compare the EFT expansion with the standard calculation of 
the fixed-width Born cross section. 

To see the correspondence with the EFT calculation, we parameterize the W mo- 
mentum as k^ = Mwv^ 1 + r M , where is a potential residual momentum (r ~ MwS, 



5 Thc longitudinal part of the propagator is cancelled by the transverse projector from the decay into 
massless fermions. 
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r ~ MwS 1 ^ 2 ), and expand P(k) in 5, including an expansion of the self-energy around 
M^y and in the number of loops, 

Ii^(k 2 ) = M^^5 n n (m ' n) , (26) 

m,n 

with 5 = (k 2 — Mw)/Mw and m denoting the loop order. The result is 



where, to make the notation simpler, we included the QCD correction A^ 3 / 2 ) from (21) 
into AW = A« + A( 3 / 2 ) instead of expanding it out, and A( 2 ) = M w (Il^ + U^Il (1 ^). 
Next we eliminate r from the numerator in (27) by completing the square and obtain 

P(r) = ! (i + n<«» - W±g 



4M 2 (V„-^-^V 4M «' VM »" 



+ T72- • (28) 



The individual terms now have a clear interpretation in the EFT formalism. The first 
term in the second line corresponds to a single insertion of the NLO terms - a kinetic en- 
ergy correction and a second-order width correction - in the non-relativistic Lagrangian 
(11) into a W line. The local term, — i/(AM^), in the second line is similar to a corre- 
sponding term in single resonance production [13], where it contributes to a production- 
decay vertex at tree level. Here this term leads to potential loop integrals with only 
one or no non-relativistic W propagator, which vanish in dimensional regularization. 
Thus, we can drop this term. In the first line of (28) we recognize the non-relativistic 
W propagator (10) multiplied by a correction to the residue. The residue correction 
originates from the expansion of the field normalization factor w defined in (12), and 
from the derivative of the renormalized one-loop self-energy, II^ 1 ' 1 ), at k 2 = M^. In an 
EFT calculation these residue corrections are not associated with the propagator, but 
they enter the matching relations of the one-loop and higher-dimension production and 
decay vertices [13]. In order to compare with the "exact" Born cross section, where these 
terms are included, we keep these residue corrections here rather than in the matching 
calculation of Section 4.1. 

The real part of n^ 1 ' 1 ) depends on the W field-renormalization convention in the 
full theory. In the following we adopt the on-shell scheme for field renormalization, 
Ren^ 1,1 ) = 0, and the pole scheme for mass renormalization. Since ImWp(k 2 ) = 
—k 2 T ( ^/M w 6(k 2 ) at one-loop due to the decay into massless fermions, it follows that 
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n( 1 - 1 ) = -iV^/Mw. Furthermore, A« = M W U^ = -irg? and A( 2 ) = M W (U^ + 
n^'^n^ 1 ' ^) = — iTyy in the pole mass renormahzation scheme, which implies RelT 2 ' ) = 
(T$/M w ) 2 , Imn( 2 '°) = -r$/M w for the renormalized two-loop self-energy at k 2 = 
M w . The QCD correction A^ 3 / 2 ) = — iT w ^ can be included into —ir$ as before. 

To compare with the "exact" Born cross section, we write (25) in this renormahzation 
scheme in the form 

k 2 - m w - r$ 2 - iM w (k 2 r$/M w + rg>) / 5 x 

p(]A — i \ L |_ o ( ) (29) 

(k 2 — m w — rf) 2 + Ml (k^/M w + r«) 2 Vm 2 ,; • 

The fixed- width prescription corresponds to replacing k 2 T w ^ /M w by T w ^ in the denom- 
inator, but not in the numerator, where the factor of k 2 arises from the integration over 

the two-particle phase space of the W decay products. In addition one drops the T w ' 
terms (since they come from Rel^ 2,0 )) and r$. Repeating the derivation of (28) with 
this modified expression we obtain 

r i r (0)2 

P(A;) fi xcd-width = Eq. (28) with A( 2 ) = -ir# - + f . (30) 

L J (k 2 - M 2 w f + M W T^ 

The additional term is purely real and does not contribute to the cut propagator ImP(fc) 
relevant to the cross-section calculation. We therefore arrive at the interesting conclusion 
that the fixed-width prescription coincides with the EFT approximation in the potential 
region up to the next-to-leading order, if Mw is the pole mass, up to a trivial term 
related to the one-loop correction T w ^ to the pole scheme decay width. 

In the calculation of the NLO correction to the forward-scattering amplitude in the 
potential region, we use the expansion (27) in (24), and drop all terms beyond NLO. This 
already accounts for all NLO corrections from the effective Lagrangian, and for some 
corrections from higher- dimension production operators with tree-level short-distance 
coefficients. Further corrections of this type come from the expansion of the squared 
matrix element &(E,r). The square of the production amplitude of two off-shell W 
bosons depends on four kinematic invariants, which we may choose to be r 2 , pi ■ r, 
k\ — M w , and k\ — M w . This choice is convenient, since all four invariants are small 
with respect to M w in the potential region. In the expansion of $(E,r) to NLO, we 
may further approximate r 2 by — r 2 , since r ~ r 2 /Mw <C \r\ and exploit that P{k\^) 
does not depend on the direction of r. We find, for the e^e\ and e~^e\ helicity initial 
states (the LL and RR combinations vanish), 



2„2 
ew 



, 11 o,2/ \ 38^ A f 2 
My + ^( S ) + Y ^))^r 



+ 0(5 2 ), 



<S> RL (E,r) = -l28n 2 a 2 ewX 2 (s){^ + 0(5 2 ). (31) 



14 



The functions 



sis 



Mf) 



X(s) = 



6M^M 



2 „2 



Mf) 



(32) 



originate from the s-channel photon and Z boson propagators. The NLO terms pro- 
portional to r 2 can be identified with tree- level production operator products O^Oljp 
and O^^Op 1 ^ 2 ^ as discussed in Section 2.3. In such calculations £(s) and x( s ) would be 
evaluated at s = 4M^. Here we keep the exact s-dependence, since this can be done at 
no calculational cost. 

Note that the coefficient functions of production operators in the EFT are determined 
by on-shell matching, which implies an expansion of amplitudes around the complex pole 
position s = + M W A rather than [30,31]. The difference cannot be neglected 
in NLO calculations. In principle the expansions (31) could have yielded terms such as 
k\ — Myy, which should be written as k\ — s + MvkA. The difference k\ — s cancels 
a resonant propagator (possibly giving rise to a production-decay operator matching 
coefficient), while the remaining My/A term must be combined with other contributions 
to the loop correction to the leading-order production vertex. This complication can be 
ignored here, since the expansion of <&(E,r) is independent of k\ 2 — up to NLO. 

The NLO correction from the potential region is now obtained by inserting the ex- 
pansions (27), (31) into (24) and performing the loop integral. The integral has an odd 
power-divergence which is finite in dimensional regularization. The LO cross section has 
already been given in (17). The NLO terms are 



r (1) 

LR, Born 



47TO! 2 

274^ 



H + 2£ 2 ( S ) + ^(s) 11m 



3/2 



Mi 



w 



+ Im 




a {1) 

w RL, Born 



8tt« 2 
27^ 



X 2 {s) Im 



3/2 



M 



w 



(33) 



Since E/M w ~ T^/M w ~ 5 and r$/M w ~ 5 every term is suppressed by 5 relative 
to the leading order as it should be. The unpolarized cross section is one fourth the sum 
of the LR, RL contributions. The factor 1/27 comes from the tree- level branching ratio 
for the final state piTv^ ud in the conversion from the forward-scattering amplitude to the 
partial cross section. As discussed above, when we use this expression to compare with 
the standard Born cross section in the fixed-width scheme, we set to zero. When we 
use the expression (33) in the complete NLO calculation including radiative corrections, 
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(1) (2) (3) 



Figure 4: Example of a two-loop diagram with one hard and one potential loop. Cut (2) 
is part of the Born cross section, but subleading as discussed in the text. 

we have to keep in mind that multiplying all terms by the product 1/27 of leading- 
order branching fractions as in (33) is actually not correct. The required modification is 
discussed in Section 3.2. 

In addition to the 5-suppressed terms from the potential region of the one-loop dia- 
grams shown in Figure 3, there is another NLO contribution from the leading terms of 
two-loop diagrams with one hard and one potential loop, which may also be associated 
with the Born cross section. An example is displayed in Figure 4. Cut (1) does not 
correspond to a four-fermion final state and must be dropped. Cut (3) corresponds to 
the interference of a tree-level production operator with the real part of a hard one-loop 
correction to a production operator. Since the s-channel diagrams do not contribute to 
the leading-power production operator, this cut is beyond NLO. Cut (2) is a contribution 
to what is usually termed the "Born cross section" corresponding to the interference of 
single and double resonant diagrams in the kinematic region where both fermion pairs 
have invariant mass of order M^. The contribution from this cut is contained in the 
imaginary parts of the hard one-loop correction to the production operators. The thresh- 
old suppression of the s-channel diagrams applies here as well, hence this contribution 
is also not relevant at NLO. 

3.2 Decay-width correction for the flavour-specific cross section 

As already noted, the expression (33) has to be modified in order to take the radiative 
correction to the decay correctly into account. In this subsection we derive the required 
modification of the formula, but note that it will not be needed for the comparison to 
the Born cross section, where radiative corrections are excluded. 

To include the loop corrections to W decay for the flavour-specific four-fermion final 
state \i~Vp, ud we have to identify contributions to the forward-scattering amplitude from 
cut two-loop W self-energy insertions and include only the appropriate cuts containing 
a muon and muon-antineutrino or up and anti-down quarks and, possibly, a photon. 
Repeating the expansion in the potential region performed in Section 3.1 for the cut dia- 
gram with flavour-specific cuts selected, one finds that in the pole mass renormalization 
and on-shell field renormalization scheme adopted here all terms in the expansion are 
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correctly treated by multiplying the totally inclusive result by the ratio of leading-order 
partial branching fractions, r^°2 P r^/[r^ ] 2 = 1/27, except for one term involving the 

insertion of A^ = — iV$ . In (33) this insertion results in part of the term involving 
r$, and is also multiplied by 1/27. We therefore have to modify this term to include 
the flavour-specific cuts correctly. At NLO we have to consider diagrams where iA^ /2 
is inserted in only one of the two H^-lines. Cutting this line produces a contribution to 
the imaginary part of the forward-scattering amplitude of the form 



Im 



% iA^ i 
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= — Im 
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where rj is the inverse propagator of the non-relativistic W boson. The first two terms 
correspond to cutting the W line to the left and right of the A^ 2 ) insertion. The flavour- 
specific final states are extracted from these cuts as discussed below (16). This amounts 
to multiplying the NLO correction (33) by the leading-order branching ratios, so these 
two terms are treated correctly by the factor 1/27. The last term corresponds to a cut two 
loop self-energy insertion, where only the cuts leading to the desired final state must be 



i(l,eui) 



taken into account. Therefore here — ImA^ = Y\J has to be replaced by T^l p ^ = T 

and T^j = F^f v) + 1.409 ^T^, respectively, to obtain the NLO cross section for the 
four-fermion final state. To implement these replacements, note that the contribution of 
the last term in (34) to the forward-scattering amplitude is of the form /F^) Im^ -*. 
We can therefore compensate the incorrect treatment of the flavour-specific cross section 
in (33) by subtracting this contribution for each W line and adding the flavour-specific 
corrections. Multiplying by the leading-order branching fraction for the second W line 
one obtains the additional NLO correction to the cross section, 



,(0) 



L 



(i) 

ud 



<0)_ 
ud 



r (l)' 

2 w 



w 



a 



(0) 



(35) 



At NLO this correction is equivalent to multiplying the imaginary part of the leading- 
order (or even next-to-leading order) forward-scattering amplitude by the one-loop cor- 
rected branching ratios T^ I l ' ) r^j L0 ' ) /[r[^ L0 ^ 2 rather than by 1/27, where r^ LO ' ) = 



,(0) 



+ rS>. 



The NLO partial decay rates are calculated in Section 4.1. 



3.3 Expansion in the hard region 

We now consider the hard contributions, which determine the matching coefficients of 
four-electron production-decay operators. As already discussed in Section 2.2, the one- 
loop diagrams shown in Figure 3 do not provide imaginary parts of the forward-scattering 
amplitude. The leading hard contributions originate from the two-loop diagrams in 
Figure 5. These diagrams are to be calculated in standard perturbation theory with no 
width added to the W propagator, but expanded near threshold. The result must be of 
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Figure 5: Two-loop cut diagrams. Symmetric diagrams are not shown. 

order a 3 , which results in a N^LO correction relative to the leading-order cross section. 
Higher-order terms in the hard region come from higher-order terms in the expansion (in 
E = y/s — 2M\y) near threshold and from diagrams with more hard loops, all of which 
are N 3 / 2 LO and smaller. 

In the hard region it is simpler to calculate the four-fermion cross section directly as 
the sum over the relevant cuts of the forward-scattering amplitude as shown in Figure 5. 
Note that this includes cutting W lines as well as diagrams with self-energy insertions into 
the W propagator. This can be interpreted as an expansion of the resummed propagator 
in the distribution sense [32,33], such as 

Mw^W xn 2 ,,2\ , T3"\ t M W T W ( 5 2 \ 

FWT^w = * Hk ~ Mw) + pv W=»Q? +0 \MJ- m 

"PV" denoting the principal value. The left-hand side arises from cutting fermion-loop 
insertions into the W propagator, but not the W lines itself. But the leading term in 
the expansion of this expression, equivalent to the narrow-width approximation, looks 
as if a W line with no self-energy insertions is cut. 

The principal-value prescription is redundant at N 1 / 2 LO, where the singularity in the 
integrand is located at one of the integration limits, and is regularized by dimensional 
regularization, which has to be supplied in any case to regulate infrared divergences 
that arise as a consequence of factorizing hard and potential regions in the threshold 
expansion. As in the potential region, the integrals are actually analytically continued 
to finite values, since the divergences are odd power divergences. The result of the 
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calculation can be written as 
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(37) 



Here the first sum extends over the diagrams as labelled in Figure 5, the second over 
the fermions / e u,d, fi, in the internal fermion loops. The explicit values of the 
coefficients arising from the diagrams hl-h3 are 

K hl = -2.35493 , K h2 = 3.86286 , K h3 = 1.88122. (38) 

The three coefficients contain the contribution of the diagrams hl-h3 shown in Figure 5 
and of the symmetric diagrams with self-energy insertions on the lower W line. K h2 
contains also the contribution of the complex conjugate of h2. The explicit expressions 
of coefficients K{ and Cf h , with h = LR, RL, for the diagrams h4-h7 are given in 

Appendix A. Similar to (32) the s-dependence of the C{ h arises trivially from photon 
and Z propagators, and we could put s = AM^ at N 1//2 LO. Since all other terms in 
(37) are energy- independent, we conclude that the leading hard contribution results in 
a constant N 1 / 2 LO shift of the cross section. 

This contribution can be interpreted as arising from a final state where one fermion 
pair originates from a nearly on-shell W decay, while the other is produced non-resonant- 
ly, either from a highly virtual W, or as in the truly single-resonant diagrams h4-h7. 
Numerical investigation reveals that the contribution from h4-h7 is rather small, below 
0.5% of the full tree cross section in the energy range y/s = 155 GeV and 180 GeV. Below 
155 GeV it becomes negative and its magnitude grows to 4% at 150 GeV. The smallness 
of the single-resonant contributions is in part due to large cancellations between the 
diagrams h4 and h5. 

The comparison with the Born cross section performed below shows that the region of 
validity of the EFT expansion is significantly enlarged, if the energy-dependent N 3 / 2 LO 
terms are included. These can only arise from the next-to-leading order terms of the 
expansion in the hard region (the expansion in the potential region produces only integer- 
power corrections in S). The energy-dependent terms are related to the next order in 
the threshold expansion of the cut diagrams in Figure 5. The computation for the 
numerically dominant diagrams hl-h3 gives 

[K a hl +K a h2 as)+K a h3 e(s)}, 

KsX 2 (s), (39) 
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where 

K a hx = -5.87912 , K a h2 = -19.15095 , K a hz = -6.18662. (40) 

Other N 3 / 2 LO corrections related to the Born cross section arise from cut three-loop dia- 
grams of the type hl-h3, but with two self-energy insertions, and of type h4-h7 with one 
insertion. This N 3//2 LO term is (almost) energy-independent and can be parameterized 
by 

«-^E^w«f- («) 

The coefficients C^ h (s) are equal to the factors multiplying in (39) and we omitted the 
small contributions from h4-h7. The calculation of the numerical coefficients K\ is non- 
trivial, since it contains products of distributions. A rough estimate of these corrections is 
a hB 2 2n ~ a hvl^w /Mw ~ 0.025 <y^ 2 \ resulting in an energy-independent contribution 
to the cross section of order 2fb. The comparison below suggests that actually it is 
significantly smaller. 



3.4 Comparison to the four-fermion Born cross section 

We compare the successive EFT approximations to the four-fermion Born cross section 
in the fixed-width scheme. We discuss only the unpolarized cross section given by (<Jlr + 
<Jrl)/4:. The relevant terms are given in (17), (33), (37), and (39). The input parameters 
are taken to be 

M w = 80.403 GeV, M z = 91.188 GeV, G M = 1.16637 • 10~ 5 GeV~ 2 . (42) 
The pole mass Mw is related to the on-shell mass through the relation (valid to 0{T W )) 

^ = M - + 2^' (43 > 



w = -,-M w = ^V 51 • ( 44 ) 



where 

„ 3 a 3G,M3 
4 s 2 w 2V2ti 

We use the fine-structure constant in the scheme, a = \f2G fJi M w s1 /re, and the on- 
shell Weinberg angle c w = M w /M z . Inserting (44) into (43), and solving the equation 
for M w , we get the following pole parameters: 

M w = 80.377 GeV, T w = 2.04483 GeV. (45) 

The value of the W width used here is the leading-order decay width (19), excluding 
the one-loop QCD correction. This is appropriate for a tree-level calculation and en- 
sures that the branching ratios add up to one. Correspondingly we set A^ 2 ) = in 
the effective-theory calculation. In Figure 6 we plot the numerical result obtained with 
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Figure 6: Successive EFT approximations: LO (long-dashed/blue), N 1,/2 LO (dash- 
dotted/red) and NLO (short-dashed/green). The solid/black curve is the full Born result 
computed with Whizard/CompHep. The N 3//2 LO EFT approximation is indistinguish- 
able from the full Born result on the scale of this plot. 

Whizard [17] for the tree- level cross-section, and the successive effective-theory approxi- 
mations. We used the fixed-width scheme in Whizard and checked that the results from 
the O'Mega [34], CompHep [18] and MadGraph [20] matrix elements agree within the 
numerical error of the Monte-Carlo integration. The large constant shift of about 100 fb 
by the N 1//2 LO correction from the hard region is clearly visible, but the NLO approx- 
imation is already close to the full Born calculation. In Table 1 we perform a more 
detailed numerical comparison, now including also the N 3//2 LO approximation. (The 
missing energy-independent N 3//2 LO terms are set to zero.) We observe that the conver- 
gence of the expansion is very good close to the threshold at yfs ~ 161 GeV, as should be 
expected. The accuracy of the approximation degrades as one moves away from thresh- 
old, particularly below threshold, where the doubly-resonant potential configurations are 
kinematically suppressed. If one aims at 0.5% accuracy of the cross section, the NLO 
approximation suffices only in a rather narrow region around threshold. Including the 
N 3//2 LO term from the first correction in the expansion in the hard region leads to a 
clear improvement both above (~ 0.1% at 170 GeV) and below threshold (~ 10% at 
155 GeV). The energy region where the target accuracy is met now covers the region of 
interest for the W mass determination (see Section 6.4). 

4 Radiative corrections 

In this section we calculate the NLO contributions that correspond to genuine loop 
corrections to four-fermion production. As outlined in Section 2.3 there are several such 
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EFT(LO) 


EFT(v / NLO) 


EFT(NLO) 


EFT(N^LO) 


exact Born 


155 


101.61 


1.62 


43.28 


31.30 


34.43(1) 


158 


135.43 


39.23 


67.78 


62.50 


63.39(2) 


161 


240.85 


148.44 


160.45 


160.89 


160.62(6) 


164 


406.8 


318.1 


313.5 


318.8 


318.3(1) 


167 


527.8 


442.7 


420.4 


429.7 


428.6(2) 


170 


615.5 


533.9 


492.9 


505.4 


505.1(2) 



Table 1: Comparison of the numerical computation of the full Born result with Whizard 
with successive effective-theory approximations. 

contributions: an electroweak correction to the matching coefficient of the leading W 
pair-production operator and to W decay; a correction from potential photons associated 
with the Coulomb force between the slowly moving W bosons; and soft and collinear 
photon effects. 

4.1 Hard corrections to production and decay 

The two hard electroweak corrections required for a NLO calculation are the one-loop 
corrections C^\ R and C^ RL in the production operator (23) and the two-loop electroweak 
W self-energy A( 2 ) see (21). We reiterate that these are conventional perturbative cal- 
culations performed in a strict expansion in a ew . In particular, in the 't Hooft-Feynman 
gauge, the propagators of the massive gauge bosons are simply given by —ig^ u /(k 2 — M 2 ) 
and the self-energy insertions are taken into account perturbatively. All fermions except 
for the top quark are treated as massless. 

Before addressing these two calculations separately, we briefly discuss the renormal- 
ization conventions for the parameters and fields of the electroweak standard model 
(SM). For a scattering amplitude, whose tree-level expression is proportional to g™ w = 
(Aira ew ) n / 2 = (Ana/ s 2 w ) n l 2 the one-loop counterterm is given by 

[tree] ( -n ^ + n 5Z e + \ ^ 5Z cxt ) , (46) 

V Sw cxt / 

where the sum extends over all external lines. As specified in (42) the three independent 
parameters of the electroweak SM are taken to be the W and Z boson mass, and the 
Fermi constant (including the electromagnetic correction to muon decay in the Fermi 
theory), while c w = Mw/Mz and a = a ew s 2 w = \plG ^M"^ s 2 w j % are derived quantities. 
Similar to the a(Mz) scheme, the (-^-scheme for defining the electromagnetic coupling 
has the advantage that the light-fermion masses can be set to zero [35,36]. The coun- 
terterm for s w is related to the W- and Z-boson self-energies. In the scheme we 
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have 

Jj^L + S Z e = 1 U t Z (0) + n^(0)-Ren^(M^) _ 6r ^ 

s w s w c w 

where 11^ is the transverse self-energy of the W boson 6 and 



a / 7 -4s 



w 
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appears in the explicit expression for the electroweak correction to muon decay, Ar (see 
e.g. [36]). For the field-renormalization counterterms SZ cxt for the external lines we use 
the conventional on-shell scheme for wave-function renormalization [36] in accordance 
with the choice made in Section 3.1 for the renormalized W propagator. In particular, 
for the VF-boson and fermion wave-function renormalization we have 

6Z w = Re^^\ p2=M2w , 5Z f = ReUf(0), (49) 

where 11^ denotes the self energy of the fermion. (Note that Re 11^(0) = 11^(0).) The on- 
shell field renormalization of the fermions ensures that no further finite renormalization 
is needed in calculating the scattering amplitude. On the other hand, since we never 
consider a physical process with external W bosons, the renormalization factor for the 
W field is purely conventional, and our final result is independent of the convention for 
5Z W . However, the matching coefficient of the production operator calculated below 
does depend on this convention. The dependence is cancelled by the dependence of (28) 
on IT 1 ' 1 ), the on-shell derivative of the renormalized one-loop self-energy, whose value 
depends on 5Zw- 



4.1.1 Production vertices 

The general method on how to obtain the matching equations needed to determine the 
short-distance coefficients of production operators has been discussed in [13]. For C^\ R 

and C ( p 

1 rl we compute the e L//ij ,e^ L — ■> W W + scattering amplitude at leading order 
in the non-relativistic approximation using dimensional regularization in d = 4 — 2e 
dimensions. This is compared to the amplitude obtained with the tree-level operator 
in the effective theory and the matching coefficient is determined to make the results 
agree. The matching coefficients thus determined are gauge invariant by construction 
provided the scattering amplitude is calculated with the external W boson momenta at 
the complex pole position. The matching prescription also includes an additional factor 
\ / 2M w w^ 1 ! 2 [13], as given in (12), for each external f2 field. However, here we depart 
from the "correct" matching procedure and omit the factor zu^ 1 / 2 , since it was already 
included in Section 3.1 (see discussion after (28)). 

6 In the conventions used here and in [13] the sum of the amputated 1PI graphs is given by (—in) 
which is the opposite sign compared to [36]. 
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Figure 7: Sample diagrams contributing to the matching of the production operator O p 
at one loop. 



The diagrams for the e~ (pi)e + (p2) — > W~(ki)W + (k2) scattering process are gener- 
ated with FeynArts [37] and the algebra is performed with FeynCalc [38]. At one loop, 
there are 65 two-point diagrams, 84 three-point diagrams and 31 four-point diagrams 
(generically counting up- type quarks, down type quarks, leptons and neutrinos), some of 
which are shown in Figure 7. Due to the simplified kinematics, many of these diagrams 
do not contribute. In fact, since the one- loop contributions are already suppressed by 
Oi e w ~ 5 it is sufficient to take the leading order in the non-relativistic expansion of the 
one- loop diagrams and to set kf, k\ to rather than to the complex pole position. 
Thus, for the W momenta we use k± — k 2 — M w v whereas the incoming lepton momenta 
can be parametrized as p\ = (Mw,p) and p 2 = (Mw, —p) with \p \ = Mw- This results 
in two simplifications. First, many diagrams vanish consistent with the fact that the 
tree-level s-channel diagrams do not contribute at leading order in the non-relativistic 
expansion. Second, the number of scales present in the loop integrals is reduced. Due to 
the simplified kinematics, all box integrals can be reduced to triangle diagrams and the 
one-loop correction to the amplitude for the process e^e^ — > W~W + takes the simple 
form 
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(50) 



expected from (23), with e^ A denoting the polarization vectors of the W bosons. (For 

h = RL, the fermion helicities are reversed.) The scalar coefficients C^ h can be obtained 
by projections of the full amplitude. Thus, we are left with the calculation of a scalar 
quantity and standard techniques for the reduction of tensor and scalar integrals can be 
applied. 

In the computation of C^ RL all poles cancel and we are left with a finite result. This 
is to be expected, since the corresponding Born term vanishes, as indicated in (13). For 
C^l R , the matching coefficient of the operator that does not vanish at tree level, the 
poles do not cancel. After adding the counterterm (46) with n — 2, it takes the form 
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(51) 



where the finite part c^lr together with the expression for C^ RL is given explicitly 
in Appendix B. For the final expression of the matching coefficient, the poles have to 
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be subtracted. However, we leave them explicit in order to demonstrate their cancella- 
tion against (double) poles from the soft contribution and poles related to initial-state 
collinear singularities. Numerically, 

= -10.076 + 0.205i (52) 

for M w = 80.377 GeV, M z = 91.188 GeV, top-quark mass m t = 174.2 GeV and Higgs 
mass M H = 115 GeV. 

The matching coefficients C^\ R and C^ RL both have a non-vanishing imaginary part. 
Taken at face value, this imaginary part contributes to the imaginary part of the forward 
scattering amplitude A and, therefore, to the total cross section. Denoting by A^ c the 
NLO contribution to A resulting from we have 

ImA^c = Im (2C^A {0) ) = 2ReC^ lmA {0) + 2ImCW ReA {0) . (53) 

However, the second term in (53) is induced by cuts that do not correspond to the final 
state we are interested in, such as the Z7 intermediate state in the fourth diagram of 
Figure 7. In fact, at leading order in the non-relativistic expansion, none of the diagrams 
that contribute to the hard matching coefficients contains either a quark or a muon. To 
obtain the flavour-specific cross section we are concerned with, we therefore have to drop 
the second term in (53) and in what follows it is always understood that we take the real 
part of the matching coefficients C^ LR and C^ RL . Recalling the discussion of cut (2) at 
the end of Section 3.1, we note that beyond NLO the situation is more complicated, as 
some of the cuts contributing to the imaginary part of the matching coefficient C p do 
correspond to the flavour-specific cross section we are interested in. 

The contribution to the cross section resulting from the NLO correction to the pro- 
duction operators is obtained by multiplying the imaginary part of A^ c by the leading 
order branching ratios. The correction to the cross section for the e^e^ polarization is 
therefore given by 

A4 1 a ) rd = ^2Re^Im4t (54) 

Because there is no interference of the helicities e R e~^ and e^e R , the coefficient C^ RL 
does not contribute at NLO. Introducing the abbreviations 

?2 p(0) -2 r (0) 

' 2M W 2 l+ 2M W 2 v ; 

for the non-relativistic propagators in the leading-order diagram, Figure 2, and jl 2 = 

/i 2 e 7E /(47r), we can rewrite (54) as 

^ = ^W f ddr l ' e 



(56) 



hard """ 27M 2 v s \ y ^ J (27r) rf ^ + 



x 2Re — 
2ir 



1 3 \ / 4M 2 



2 2eJ\ ^ I "- LH 



W \ + c (l,fin) 



The unintegrated form of the result is given to make the cancellation of the e-poles 
against other contributions computed in the following subsections more transparent. 
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Figure 8: Diagrams contributing to the virtual correction C^p h at one loop. 



4.1.2 Decay corrections 

Next we discuss the electroweak correction to the matching coefficient A. In the pole 
mass and on-shell field renormalization scheme A^ 2 ' ew ^ = — iT^ 1 '^ = iM w Imll^ 2 ' ). The 
cuts of the 2-loop electroweak W self-energy consist of two parts, corresponding to the 
virtual and real hard corrections to the W pole decay width. Dealing with the total 
cross section, we only need the sum of these two. However, we also have to discuss how 



to obtain results for the flavour-specific process e 



/i v^udX. To aid this, we will 



discuss the virtual and real corrections separately, starting with the former. 

The virtual one-loop correction to the pole-scheme decay width into a single lepton 
(I) or quark (h) doublet can be written as 
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where the tree-level widths in d dimensions are V 
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(57) 

= a ew M w /12 + 0(e) and 

r^°j = 3T^. The calculation of cfl involves the evaluation of the diagrams 
depicted in Figure 8 with obvious modifications for the leptonic decay. After adding the 
counterterm (46) with n = 1 we obtain 
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where for the leptonic (hadronic) decay we have to set the electric charges to Qf — 
— 1,0/ — (Qf — 2/3, Qf = —1/3 ). The finite parts c^'^ of the matching coefficients 
are given explicitly in Appendix B. Numerically, 
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= -2.709 - 0.552 i, 



„(i,fin) 

-d,h 



= -2.034 - 0.597 i, 



(59) 



for M w = 80.377 GeV, M z = 91.188 GeV, m t = 174.2 GeV, and M H = 115 GeV. 

To this we have to add the correction due to hard real radiation of a single photon. 
Since the corresponding soft corrections vanish, the hard real corrections are equivalent 
to the real corrections evaluated in the standard electroweak theory and their calculation 
is straightforward. We compute the bremsstrahlung diagrams and integrate the squared 
amplitude (divided by 2M W ) over the <i-dimensional phase-space [39]. The expression 
thus obtained contains infrared (double) poles which cancel the poles in (57) and we 
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Figure 9: First and second order Coulomb correction. 



are left with finite expressions for the flavour-specific leptonic and hadronic matching 
coefficients. Including the (two-loop) QCD correction to the hadronic decay, they read 
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(60) 



Strictly speaking, for the computation of these matching coefficients we have to expand 
around the complex pole s and not around M^. However, the difference in the width is 
of order a 3 and thus beyond NLO [2]. 



4.2 Coulomb corrections 

The exchange of potential photons with energy /c ~ MyyS and three-momentum k ~ 
M w y/5, shown in Figure 9, corresponds to insertions of the non-local four-boson inter- 
actions in the effective Lagrangian (11). These insertions can be summed to all orders 
in terms of the Green function G c (f, f E) of the Schrodinger operator — V 2 /Mw — ot/r 
evaluated at r = r' = 0. Using the representation of the Green function given in [40], 
we obtain [41] 




iA nnih>m u= - unai w al^[- EJ [f^ )+Al = 1 }. (CI) 
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where ip(x) is Euler's psi-function, and a subtraction-scheme dependent real constant 
that drops out in the cross section has been omitted. The diagram with no photon 
exchange is not included in this expression, since it corresponds to the leading-order 
amplitude (15). The logarithm constitutes a a/y/5 ~ V5 correction relative to the 
leading-order scattering amplitude (15). The expansion of the psi-function in a results 
in an expansion in powers of V5. Thus, the Coulomb correction up to NLO reads 
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(62) 
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(im) (mml) (mm2) 



Figure 10: Soft-photon diagrams in the effective theory: Initial-initial state interfer- 
ence (ii), initial-intermediate state interference (im) and intermediate- intermediate state 
interference (mm). Symmetric diagrams are not shown. 



This contributes only to the LR helicity cross section, since the production operator at 
the vertices in Figure 9 is the leading order one (13). Directly at threshold (E = 0) the 
one-photon exchange N 1 ^ 2 LO term (the logarithm in (62)) is of order 5% relative to the 
leading order. Two-photon exchange is only a few-permille correction, confirming the 
expectation that Coulomb exchanges do not have to be summed to all orders due to the 
large width of the W boson. The one and two Coulomb-exchange terms have already 
been discussed in [10,11]. 

4.3 Soft-photon corrections 

We now turn to the radiative correction originating from soft-photon exchange. These 
are 0(a) contributions to the forward-scattering amplitude, and correspond to two- 
loop diagrams in the effective theory containing a photon with momentum components 
Qo ~ \q\ ~ MyyS- The relevant Feynman rules are given by the coupling of the soft 
photon to the Q± fields in the PNRQED Lagrangian (11) and to the collinear electrons 
and positrons contained in the SCET Lagrangian. The latter is simply the eikonal 
coupling ±ien M , where n M is the direction of the four- momentum of the electron or 
positron. The topologies contributing to the two-loop forward-scattering amplitude are 
shown in Figure 10. The W^-boson vertices are leading-order production vertices, hence 
at NLO the soft correction applies only to the left-right e~e + helicity forward-scattering 
amplitude. Note that (mm2) is not a double-counting of the Coulomb-exchange diagram 
in Figure 9, since the two diagrams refer to different loop momentum regions. 

It is well known that for the process e~e + — > W~W + — > /1/2/3/4 the soft-photon 
corrections related to the final state cancel for the inclusive cross section [42,43]. The 
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diagrams of type (im) in Figure 10 cancel pairwise when the sum over incoming positrons 
and electrons is performed. The sum of the diagrams of the form of (mml) and (mm2) 
cancels after the loop integrals are performed. Therefore the sum of all diagrams where 
a soft photon couples to an Q line vanishes. In the effective theory this cancellation can 
be seen from the outset, since it follows from the particular form of the leading coupling 
of a soft photon to non-relativistic W bosons in the effective Lagrangian (11), which 
involves only A° s (t,0). Since the residual gauge invariance of the effective Lagrangian 
allows one to set the time-like component of the photon field to zero, at leading order 
the couplings can be removed from the Lagrangian. 

Therefore the soft-photon correction in the effective theory is given by the initial- 
initial state interference diagrams. However, diagram (ii2) leads to a scaleless integral 
which vanishes in dimensional regularization, and diagram (ii3) and the symmetric dia- 
gram are proportional to pi ~ and p\ ~ 0, respectively. The only non-zero diagram is 
(iil) and the corresponding crossed diagram. The sum of the two diagrams evaluates to 
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The double e-pole in (63) cancels against the pole in the hard matching coefficient; the 
single pole can be factorized into the initial-state electron (positron) structure function 
as shown in Section 5. Subtracting the pole part of the integrand (63) before performing 
the integration, one obtains 
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(64) 



As before, the r° integration has been performed by closing the r° integration contour in 
the upper half-plane and picking up the pole at r° = E — f 2 / (2Myy) +iF^) /2. Because of 
the absence of soft corrections related to the final state, at NLO the soft corrections to the 
flavour-specific process (1) can be obtained by multiplying the soft two-loop contributions 
to the forward-scattering amplitude by the leading-order branching ratios, thus 



Aa W - — ImAi (1) 



(65) 



As a check, we also calculated the soft corrections directly for the process (1) and found 
agreement with the simpler calculation of the forward-scattering amplitude. 



29 




Figure 11: Collinear-photon diagrams in the effective theory. Two symmetric diagrams 
are not shown. 



4.4 Collinear-photon corrections 

Finally we consider collinear-photon corrections, corresponding to photon energies of 
order Myy, and photon virtuality of order My^Ty/. The four-momentum of the photon 
is proportional to the initial-state electron or positron momentum. The collinear photon 
couplings arise from the SCET Lagrangian, while their couplings to the W bosons is 
encoded in the collinear Wilson lines in the production operators. The diagrams corre- 
sponding to NLO contributions are shown in Figure 11. As discussed in [13] all these 
diagrams are scaleless for on-shell, massless initial-state particles. However, we shall 
have to say more about collinear effects in Section 5, when we include the resummation 
of large initial-state radiation logarithms. 



4.5 Summary of radiative corrections 

The radiative correction to the next-to-leading order cross section is given by the sum 
of the corrections (56), (65), (62), (35) computed in the previous sections, 
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Recall that this refers to the e^e^ helicity initial state, while there are no radiative 
corrections to the other helicity combinations at NLO. The radiative correction to the 
unpolarized cross section is one fourth of the LR contribution. 

Because of the approximation m e = 0, the cross section is not infrared-safe, as can 
be seen by summing the four contributions. The Coulomb and decay corrections are free 
of infrared singularities. For the sum of the soft (63) and hard (56) terms we obtain the 
following expression: 
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The cross section a L ' R is a "partonic" cross section. It should be convoluted with the 
electron (positron) distribution function, which contains the infrared effects associated 
with the electron mass scale. In the following section we discuss how the partonic cross 
section is transformed to the infrared-finite physical cross section. 

5 Initial-state radiation 

The remaining e-poles in (67) are associated with emission of photons collinear to the in- 
coming electron or positron, and can be factorized into the electron distribution function 
rf e s , in terms of which the physical cross section a reads [44,45] 

a h (s) = f d Xl f dx 2 r™( Xl )r™(x 2 ) ar( Xl x 2 s) . (68) 
Jo Jo 

Here <3"^ s (s) = c^,Bom(s) +^| 1 1 ms( s ) * s our resu lt f° r the NLO helicity-specific cross section 
after adding the Born cross section from Section 3 and the radiative correction from (66) 
with the infrared e-poles minimally subtracted. The partonic cross section depends on 
the scales Q = {M W) E,T W } and the factorization scale ji. The electron distribution 
function in the MS scheme depends on /i and the very- long distance scale m e . The 
physical cross section is independent of /j, and includes the electron-mass dependence 
up to effects suppressed by powers of m e /Q. By evolving the electron distribution from 
the scale m e to the scale Q, one sums large collinear logarithms a ni In™ 2 (Q 2 /ml), with 
rii = l,...,oo, n 2 = l,...,ni from initial-state radiation of photons to all orders in 
perturbation theory. A NLO calculation of the partonic cross section should go along 
with a next-to-leading logarithmic approximation, where all terms with n 2 = n\ and 
n 2 — ni — 1 are summed. Note that here we do not attempt to sum logarithms of 
M w /T w , which are less important, although the effective-theory formalism is ideally 
suited for this summation as well. 

Unfortunately the structure functions T^(x) available in the literature do not corre- 
spond to the MS scheme and sum only leading logarithms a n ln n (Q 2 /ml). To convert our 
result <x^ is (s) to this scheme and sum the leading-logarithmic initial-state radiation effects 
we proceed as follows: first, using the expansion r^ e s (x) = 5(1 — x) + T™e^ l \x) + 0(a 2 ), 
we compute the scheme-independent NLO physical cross section without summation of 
collinear logarithms, 

< LO (*) = + *Sb(8) + 2 f dxTf^\x) a h , B Uxs). (69) 

Jo 

Then, by comparing this to the corresponding equation in the conventional scheme, 

< LO (s) = <T h , Boin (s) + *£L» + 2 f dxT^\x) a h>Boin (xs), (70) 

Jo 
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we determine 0-JJ onv ( s ), and hence a c h onv (s) = <7 Morn (s) + 5"^ onv (s). Finally, we calculate 
the initial-state radiation resummed cross section 

a h (s)= [ d Xl [ dx 2 rf e (x 1 )r^(x 2 )ar v (xix 2 s) (71) 
Jo Jo 

in the conventional scheme for the electron (positron) distribution functions. Note that 
since the Born cross section for the RL helicity combination is already a NLO effect, the 
scheme conversion must be performed only for h = LR. For h = RL we simply have 
&%L(s)=a%i(s) = <T RL>Bmn (s). 

Step 1: Calculation of the fixed-order physical cross section crf^°(s). Rather than cal- 
culating the last term on the right-hand side of (69), we compute directly the radia- 
tive correction to the physical cross section, by converting ^"^^^(s), where the 
collinear divergences are regulated dimensionally, into the expression when the electron 
mass itself is used as the regulator. 

In the presence of the new scale m e <C Tw, E, My/ there are two new momentum 
regions that give non-zero contributions to the radiative corrections. They correspond 
to hard- collinear photon momentum (q° ~ M w , q 2 ~ m 2 ) and soft- collinear photons 
(q° ~ T w , q 2 ~ m 2 e Y 2 v jM V7 )? The corresponding loop integrals are scaleless when 
m e = 0; for m e ^ 0, they supply the difference 



- = Ai + Aa^ oU . (72) 



In other words cr^(s) is the sum of the four contributions in (66) plus those from the 
two new momentum regions. 

Only a small subset of all the radiative correction diagrams has hard- or soft-collinear 
contributions, namely those containing a photon line connecting to an external electron 
or positron. The topology of the soft-collinear and hard-collinear diagrams is identical to 
the (ii) and (im) diagrams in Figure 10, and to the diagrams in Figure 11, respectively. 
The calculation is straightforward. In each region we simplify the integrand by neglecting 
all small terms, since the leading-order term in the expansion in each region is sufficient. 
The soft-collinear correction is 
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7 The existence of two collinear momentum regions is related to the fact that the W pair-production 
threshold region probes the electron distribution function near x = 1, where hard-collinear real radiation 
is inhibited. 
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the hard-collinear correction 
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The structure of the logarithms makes it clear that the two contributions arise each from 
a single scale, [i ~ m e T\y / 'M\y and fi ~ m e , respectively. Adding (66), (73), (74), and 
making use of (67) results in the factorization-scheme independent radiative correction 
to the physical cross section, 
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After performing the r-integral we may set d to four and obtain a finite result. As 
expected the e-poles have cancelled, but the infrared- sensitivity of the cross section is 
reflected in the large logarithms \a(2M w /m e ). 



Step 2: Calculation of a c £^(s). Comparing the right-hand sides of (69) and (70), we 
obtain the radiative correction to the conventional "partonic" cross section 
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where Tee'^\x) is the 0(a) term in the expansion of the conventional electron structure 
function provided in [46,47]. In the notation of [47] we employ the structure function 
with j3 exp = f3 s = f3 H = f3 e = ^ (21n(y / s/m e ) — 1). To calculate the subtraction term in 
(76) it is sufficient to approximate y/s = 2M\y in the expression for /3 e , to set o~LR,Bom(xs) 
to the leading-order Born term (17) with the replacement of E by E — M w (l — x), and 



to use T^ 1 \x) in the limit x — > 1, 
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We then reintroduce the integral over r, and exchange the r- and x-integration to obtain 
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which shows that ^°r v ( s ) * s fr" ee from the large electron mass logarithms. To obtain the 
final form in (78) we have shifted the integration variable r$ to E — tq. Summing (75) 
and (78), and performing the r-integration, gives the final result for the next-to-leading 
order radiative correction to the conventional "partonic" cross section 
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S'tep 5: Computation of the resummed cross section. The summation of collinear loga- 
rithms from initial-state radiation is completed by performing the convolution (71) using 
the Born cross section and the radiative correction (79) together with the electron struc- 
ture functions from [46,47]. This constitutes our final result, which we shall discuss in 
detail in the following section. 



6 NLO four-fermion production cross section 

We now present our NLO predictions for the total cross section of the process e~e + — > 
yTVpudX and assess the theoretical error on the VT^-mass measurement due to the un- 
certainties in the cross-section calculation. 

6.1 Input parameters and summation of VK-width corrections 

In addition to the input parameters (42) used for the comparison of the tree cross section 
we use a s = af^(80.4 GeV) = 0.1199 and the masses 

m t = 174.2 GeV, M H = 115 GeV, m e = 0.51099892 MeV. (80) 

We use the fine structure constant a in the G M scheme everywhere including the initial- 
state radiation. With these input parameters we obtain from (60) the numerical value 
of the W width to NLO, 

rv = 3 (rg + rSj"°) + 2 (rg? fc + r£j">) s QCD = 2.09201 GeV. (si) 
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Note that we have chosen to multiply not only the leading order, but also the electroweak 
correction to the hadronic decay by the factor <5 QCD defined in (20). In the numerical 
results below we will resum the full NLO width (81) in the effective-theory propaga- 
tor (10), that is we do not perform an expansion of the propagator in the perturbative 
corrections to the matching coefficient A. We now describe how the formula for the NLO 
cross section must be modified to accomplish this summation of the width corrections. 
Readers not interested in this technical detail may move directly to the next subsection. 

Leaving A = — iTw unexpanded amounts to setting to zero in the NLO tree 
cross section (33) and to replacing by 1TV wherever it appears. Some care has then 
to be taken in order to obtain the correct cross section for the flavour-specific four- 
fermion final state from the calculation of the forward-scattering amplitude. Cutting the 
effective-theory propagator leads to a factor 

M f W „ ■ (82) 

fa " dfc) 2 + ~f 

analogously to (16). In the direct calculation of the four-fermion production cross sec- 
tion the numerator arises from integrating over the two-body decay phase space, which 
yields the leading-order partial width. Hence, we have to multiply all contributions 
to the forward-scattering amplitude with two cut effective-theory propagators (the po- 
tential contributions in Section 3.1, the Coulomb and soft radiative corrections, and 
the contribution from the one-loop correction to the production operator) by a factor 
F % T ul/ T w instead of the factor T^l.T^/lT^} 2 = 1/27 used in the tree level analy- 
sis. In the calculation of the matching coefficient of the four-electron production-decay 
operator performed in Section 3.3 the self-energy insertions on one of the two W lines are 
treated perturbatively, and the decay subprocess is already correctly included at lowest 
order, while the other W is effectively treated in the narrow-width approximation 



m w t w ry sn 2 n/r 2 

(P - M^f + M^T 2 W 71 T w 



2v2 lf2p2 n^Stf-M^). (83) 



To obtain the correct flavour-specific final state we therefore have to include a single 
prefactor T^__^ _ p /ry or T^ + ^ u ^/T\y, depending on the W charge. As shown in 
Table 2, with these prescriptions the N 3 / 2 LO effective-theory approximation and the 
full Born cross section (in the fixed-width definition now using (81)) are again in very 
good agreement, similar to the earlier comparison, where only was resummed in the 
propagator. 

As already mentioned the electroweak radiative corrections are correctly treated by 
multiplying the inclusive forward-scattering amplitude by T^_^^_-T^ + ^ u ^/Ty [/) except 
for the correction to W decay itself. These contributions are included by adding the 



decay correction 



'p(l,eu>) p(l,etu)' 

A*£> = I ^ + -*r ) - ,0) (<*> 



decay \ p(0) ^ p(0) 



35 





er(e e + — > /x iyud)(fb) 


y/s \GeV] 

V L J 


EFT Tree (NLO) 


EFT Tree (N 3 / 2 LO) 


exact Born 


155 


42.25 


30.54 


33.58(1) 


158 


65.99 


60.83 


61.67(2) 


161 


154.02 


154.44 


154.19(6) 


164 


298.6 


303.7 


303.0(1) 


167 


400.3 


409.3 


408.8(2) 


170 


469.4 


481.7 


481.7(2) 



Table 2: Comparison of the numerical computation of the full Born result with Whizard 
with successive effective-theory approximations as in Table 1, but now the NLO decay 
width T w as given in (81) is used. 

instead of (35). The QCD corrections up to order a 2 s are included in a similar way. 
Because of the large NLO corrections to the tree cross section and the large effect of 
ISR, it is sensible to apply the QCD decay correction to the full NLO electroweak cross 
section. This amounts to multiplying T^j, T^ w ^ by the radiative correction factor 
5qcd as given in (20), wherever they appear, which is consistent with the definition of 
the NLO W width (81). If in addition we also account (approximately) for the QCD 
decay correction to the non-resonant contributions from Section 3.3, this is equivalent 
to multiplying the entire NLO electroweak cross section by 5qcd and using the QCD 
corrected width (81) as will be done in the following analysis. 

6.2 NLO four-fermion production cross section in the effective 
theory 

The convolution of the "partonic" cross section with the electron structure functions 
contains integrations over partonic center-of-mass energies far below threshold, where 
the effective field theory approximation is not valid. The EFT calculation should be 
matched to a full cross section calculation below some cms energy, say y/s = 155 GeV, 
where for the full calculation a Born treatment is sufficient, because the cross section 
below threshold is small. Since the N 3,/2 LO EFT approximation to the Born cross sec- 
tion provides a very good approximation (except significantly below threshold), we have 
found it more convenient to replace the EFT approximation to the Born cross section 
convoluted according to (71) by the full ISR-improved Born cross section as generated 
by the Whizard program [17] rather than to perform this matching. To this we add 
the NLO radiative correction (79) (replacing the leading-order cross section by the 
full Born cross section cr Born in the decay correction (84)), which we also convolute 
with the electron distribution functions. Here we simply cut off the integration region 
\JxiX 2 s < 155 GeV. The dependence on this cut-off is negligible. Lowering it from 
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cr(e e 


f — > i/^ud X) (fb) 




x/s IGeV] 


Born 


Born(ISR) 


NLO 


NLO(ISR-tree) 


158 


61 67(2) 


4^ 64(91 


4Q 19(9) 


50 02(2) 






[-26 0%1 


[-20 2%1 

1 jZjVJ .jZj /(J 1 


[-18 9%1 

1 ± (J . U /\)\ 


161 


154 19(61 


108 60(4) 


117 81(5) 


120 00(5) 






[-29 6%1 


[-23 6%1 


[-22.2%] 


164 


303.0(1) 


219.7(1) 


234.9(1) 


236.8(1) 






[-27.5%] 


[-22.5%] 


[-21.8%] 


167 


408.8(2) 


310.2(1) 


328.2(1) 


329.1(1) 






[-24.1%] 


[-19.7%] 


[-19.5%] 


170 


481.7(2) 


378.4(2) 


398.0(2) 


398.3(2) 






[-21.4%] 


[-17.4%] 


[-17.3%] 



Table 3: Two NLO implementations of the effective-theory calculation, which differ by 
the treatment of initial-state radiation compared to the "exact" Born cross section with- 
out (second column) and with (third column) ISR improvement. The relative correction 
in brackets is given with respect to the Born cross section in the second column. 

to 155 GeV to 150 GeV (140 GeV), changes the cross section at \fs = 161 GeV from 
117.81 fb to 117.87 fb (117.91 fb), while the dependence on the cut-off for higher cms 
energy is even smaller. 

Our result for the NLO four-fermion cross section is shown in Table 3. The impact of 
radiative corrections is seen by comparing the exact Born cross section (second column, 
identical to the last column in Table 2), the ISR- improved Born cross section (third 
column) and the NLO result (fourth column). As is well-known initial-state radiation 
results in a large negative correction (about 25%). The size of the genuine radiative 
correction is best assessed by comparing the "NLO" column to the "Born(ISR)" column 
and thus seen to be about +8%. Given that we aim at a theoretical accuracy at the 
sub-percent level, this is an important effect. We shall discuss below, in Section 6.4, an 
estimate of the remaining uncertainty of the NLO cross section. 

One uncertainty is related to the fact that the conventional implementation of ISR 
sums only leading logarithms, whereas a NLO calculation of the partonic cross section 
should be accompanied by a next-to-leading logarithmic resummation. Thus rather 
than convoluting the full NLO partonic cross section with the structure functions as done 
above and indicated in (71), one could equally well convolute only the Born cross section, 
and add the radiative correction without ISR improvement, as done in some previous 
NLO calculations [6,15]. Although we favour the first option, the two implementations 
are formally equivalent, because the difference is a next-to- leading logarithmic term. 
We therefore consider this difference as an estimate of the uncertainty induced by the 
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er(e e + — > /i z/^-urf X)(fb) 


v^[GeV] 


Born 


NLO(EFT) 


ee4f [15] 


DPA [15] 


161 


150.05(6) 


104.97(6) 


105.71(7) 


103.15(7) 


170 


481.2(2) 


373.74(2) 


377.1(2) 


376.9(2) 



Table 4: Comparison of the strict electroweak NLO results (without QCD corrections 
and ISR resummation). 



missing next-to-leading logarithmic evolution of the structure functions. To assess this 
uncertainty, in the fifth column of Table 3 we show the NLO cross section based on the 
expression 

fflSR-trcc(s) = / dx l I dx 2 r^(Xi)r^(x2)<rBorn(^1^2S) + ^fonv ( S ) > ( 85 ) 

Jo Jo 

where the NLO correction to the "partonic" cross section, ffW y (s), is given in (79) (with 
1/27 replaced by T^l p T^/T^)- The comparison of the last and second-to-last columns 
of Table 3 shows that the difference between the two implementations of ISR reaches 
almost two percent at threshold and is therefore much larger than the target accuracy in 
the per-mille range. The difference between the two implementations becomes smaller 
at higher energies and is negligible at y/s = 170 GeV. The impact of this difference on 
the accuracy of the PF-mass measurement will be investigated further in Section 6.4. 

6.3 Comparison to the full four-fermion calculation 

We now compare the NLO prediction of the four-fermion production process (1) obtained 
with the effective-theory method to the full NLO calculation performed in [15] in the 
complex mass scheme. For this comparison, we adjust our input parameters to those 
of [15], 

M w = 80.425 GeV, T w = 2.0927 GeV, m t = 178GeV, a s = 0.1187, (86) 

and use a(0) = 1/137.03599911 in the relative radiative corrections as in [15]. We 
first compare the strict electroweak NLO calculation, i.e. the cross section without the 
QCD correction <5qcd and without initial-state radiation beyond the first-order term. In 
the effective-theory calculation the corresponding radiative correction is given by (75) 
omitting the second-order Coulomb correction and the factor <5qcd in the decay width. 
In Table 4 the EFT result and the result of [15] are shown in the columns labelled 
"NLO(EFT)" and "ee4f", respectively. For comparison we also show the results for the 
Born cross section and in the double-pole approximation ( "DPA" ) in the implementation 
of [6] as quoted in [15]. The main observation is that the difference between the EFT 
and the full four-fermion calculation is only 0.7% at y/s = 161 GeV and grows to about 
1% & ty/s = 170 GeV. 
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er(e e + — > /i z/^-urf X)(fb) 


v^[GeV] 


Born(ISR) 


NLO(EFT) 


ee4f [15] 


DPA [15] 


161 


107.06(4) 


117.38(4) 


118.12(8) 


115.48(7) 


170 


381.0(2) 


399.9(2) 


401.8(2) 


402.1(2) 



Table 5: Comparison of NLO results with QCD corrections and ISR resummation in- 
cluded. 



Next, in Table 5, we compare to the full result including the QCD correction and the 
resummation of ISR corrections with [15]. Here we implement the QCD correction as 
in [15] by multiplying the entire electroweak NLO result by the overall factor (1 + a s /n). 
Furthermore, we include ISR corrections only to the Born cross section as in (85), in 
agreement with the treatment of [15]. Again the second-order Coulomb correction is set 
to zero, because [15] does not include any two-loop effects. As before, the Table shows the 
two NLO calculations, the Born cross section (now ISR improved) and the double-pole 
approximation. The discrepancy between the EFT calculation and the full four-fermion 
calculation is around 0.6% at threshold. The EFT approximation is significantly better 
than the double-pole approximation directly at threshold, while at higher energies the 
quality of the DPA improves relative to the EFT approximation, since no threshold 
expansion is performed in the DPA. 

6.4 Theoretical error of the Mw determination 

The W mass will probably be determined by measuring the four-fermion production 
cross section at a few selected cms energies near the W pair-production threshold. 
In this section we estimate the error on the W mass from various sources of theo- 
retical uncertainty. To this end we assume that measurements Oi will be taken at 
y/s = 160, 161, 162, 163, 164 GeV, and at y/s = 170 GeV, and that the measured values 
coincide with our NLO calculation (labelled "NLO (EFT)" in Table 3) corresponding to 
the W pole mass Mw = 80.377 GeV. We denote by Ei(SMw) the cross section values at 
the six cms energy points for any other theoretical calculation of four-fermion production 
as a function of the input W mass 80.377 GeV + 5M W , and determine the minimum of 

i=i * 

For simplicity we assume that each point carries the same weight, so Oi = a is an 
arbitrary constant of mass dimension —2. (We checked that a more realistic assignment 
Oi ~ \f6~i does not lead to significantly different results.) The value of 5M W at which 
X 2 {SMy\r) attains its minimum provides an estimate of the difference in the measured 
value of Mw due to the different theoretical cross section inputs, Oi and Ei. For instance 
if Ei(5M w ) is the ISR-improved Born cross section (labelled "Born(ISR)" in Table 3), we 
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obtain SM W = —201 MeV, which tells us that comparing measurements to a theoretical 
calculation without the genuine radiative corrections would result in a value of M w which 
is about 200 MeV too low. The NLO calculation is therefore crucial for an accurate My/ 
determination. Next we attempt to estimate whether it is accurate enough. 

Treatment of initial-state radiation. A look at the last two columns of Table 3 re- 
veals that two different implementation of ISR, which are formally equivalent at the 
leading-logarithmic level, can lead to differences in the predicted cross section of 2% 
at \fs = 161 GeV, where the sensitivity to My/ is largest. We take this as a measure 
for the uncertainty caused by the missing next-to-leading logarithmic corrections to the 
structure function. To estimate the error on M w caused by this uncertainty we apply 
the procedure discussed above and find 

[8M W ] 1SR ~ 31 MeV. (88) 

This large error could be avoided by measuring the cross section predominantly around 
170 GeV rather than around 162 GeV, but the sensitivity to My/ is significantly smaller 
at higher energies (see Figure 12 below). Thus, this error should be eliminated by a 
consistent treatment of the electron structure functions at the next-to-leading logarithmic 
level, in which all NLL corrections are taken into account by convoluting the NLO cross 
section with the NLL structure functions. A related effect concerns the choice of scheme 
and scale of the electromagnetic coupling. The difference in the cross section between 
using cc(0) and a in the G^-scheme in the radiative correction (including, in particular, 
initial-state radiation) is about 1%, which translates into another error of about 15 MeV 
in the W mass. The scale ambiguity of the coupling used in initial-state radiation can be 
resolved only in the context of a next-to-leading logarithmic resummation which takes 
the evolution of a between m e and Yy/ into account. On the other hand, the typical 
scales in the short-distance cross section are at least rs 2 GeV, so that a in the 
scheme is more appropriate than the low-energy electromagnetic coupling in the radiative 
correction to the short-distance cross section, since it is numerically close to the running 
coupling at 2 GeV. 

Uncalculated corrections to the "partonic" cross section. The leading missing higher- 
order terms in the expansion in a and 8 are N 3 ^ 2 LO corrections to the forward-scattering 
amplitude from four-loop potential diagrams (third Coulomb correction), three- loop di- 
agrams with two potential loops and one soft loop (interference of single- Coulomb and 
soft radiative corrections), two- loop potential diagrams with 0(a) matching coefficients 
or 0(5) higher-dimensional production operators, and the 0(a) correction to the match- 
ing coefficients of the four-electron production-decay operators. The latter is expected 
to be the largest of these contributions, in particular since the non-resonant N 1,/2 LO 
contributions are large at the Born level (~ 40% at threshold, see Table 1). Presum- 
ably, this contribution is also the origin of the 0.6% difference between the EFT result 
"NLO(EFT)" and the full four-fermion calculation "ee4f" at ^ = 161 GeV in Table 5. 
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A rough estimate of this correction to the helicity-averaged cross section is 



A*~/C, (89) 

where /C is an s-independent constant of order 1. (In fact, if we attributed the difference 
between our calculation ("NLO(EFT)") and that of [15] ("ee4f") at ^ = 161 GeV 
exclusively to this contribution, we would obtain /C = 0.96.) Thus, we choose /C = 1, 
add (89) to the "NLO(EFT)" calculation, and minimize the x 2 function. From this we 
obtain an error 

[6M w ] non _ res «8MeV. (90) 

The second largest uncalculated correction to the partonic cross section is expected to 
come from diagrams with single- Coulomb exchange and a soft photon or a hard correction 
to the production vertex. A naive estimate of the sum of the two terms is 

n-W _ An {1) _ An- (1) 

A- LR Coulomb ^O decay (!) 

Act = ^ Ao-^ o ; lomb , (91) 

where the quantities involved have been defined in Section 4. Estimating the correspond- 
ing uncertainty on the W mass as before, we find 

[<5M w ] Coulombx(hard+soft) « -5MeV. (92) 

Adding the two errors we conclude that the uncertainty on Mw due to uncalculated 
higher-order effects in the effective field theory method is about 10 — 15MeV. Thus, to 
reach a total error of ~ 6 MeV requires the inclusion of at least some N 3/2 LO corrections 
in the EFT approach. The larger of the two errors estimated above, due to the elec- 
troweak correction to production-decay operator, can be removed by using the full NLO 
four-fermion calculation, where this correction is included. 



Summary. The discussion above is summarized in Figure 12, where we plot k = 
cr(s, Mw + 5Mw) / a(s, Mw) for different values of 5Mw as function of the cms energy, a 
being our NLO result, "NLO(EFT)". The relative change in the cross section is shown 
as dashed lines for 5M W = ±15, ±30, ±45 MeV. The shape of these curves shows that 
the sensitivity of the cross section to the W mass is largest around the nominal threshold 
y/s ~ 161 GeV, as expected, and rapidly decreases for larger y/s. (The loss in sensitivity 
is partially compensated by a larger cross section, implying smaller statistical errors of 
the anticipated experimental data.) 

The dark-shaded area in Figure 12 corresponds to the uncertainty on the cross section 
from (91), while the light-shaded area adds (linearly) the uncertainty from (89). The 
theoretical error decreases with y/s, since Aa in (89) is roughly energy- independent, 
while a increases. The largest current uncertainty is, however, due to ambiguities in the 
implementation of ISR. The solid (red) curve gives the ratio of the two different imple- 
mentations of ISR, NLO (EFT) vs. NLO(ISR-tree), both evaluated at M w = 80.377 GeV. 
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Figure 12: W-niass dependence of the total cross section. All the cross sections are 
normalized to a(s, M w = 80.377 GeV). See text for explanations. 

As mentioned above, we do not consider this as a fundamental problem, since this un- 
certainty can be removed with further work on a next-to-leading-logarithmic ISR resum- 
mation that will be required for many other processes at a high-energy e~e + collider as 
well. 



7 Conclusion 

We performed a dedicated study of four-fermion production near the W pair-production 
threshold in view of the importance of this process for an accurate determination of 
the P^-boson mass. Our theoretical study of radiative and finite-width corrections was 
motivated by a corresponding experimental study [3] which showed that the planned 
high-luminosity linear collider might allow a measurement of Mw with an error of only 
6 MeV from the threshold region. Our calculation, and the good agreement with the full 
NLO four-fermion cross section calculation of [15], demonstrates that accurate theoretical 
calculations are feasible and available in the threshold region. With regard to the mass 
determination, we find: 

• A resummation of next-to-leading collinear logarithms from initial-state radiation 
is mandatory to reduce the error on M w below the 30 MeV level. 
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• The NLO partonic cross-section calculation in the effective theory approach implies 
a residual error of about 10 - 15 MeV on M w . The largest missing N 3 / 2 LO effect is 
probably due to the electroweak correction to the (non-resonant) production-decay 
vertex, which is included in the full NLO four-fermion calculation, and can thus 
be eliminated. 

It is forseeable that both items can be removed, so we conclude that there is no funda- 
mental difficulty in reducing the theoretical error in the W mass determination from the 
threshold region to about 5 MeV. 

The calculation presented here is also the first NLO calculation of a realistic process 
in unstable-particle effective theory, since [12, 13] discussed the case of a single resonance 
in a gauged Yukawa model. Comparison of our results for four-fermion production 
with numerical integrations of the Born matrix elements and the radiative correction 
shows good convergence of the EFT expansion, and very good agreement once the first 
subleading term in each essential region (potential/resonant, hard/non-resonant) is in- 
cluded. The EFT approach provides a consistent treatment of finite-width effects that 
can in principle be extended systematically to higher orders. Our final results take the 
form of compact analytic formulae, which has to be compared to the numerical and 
technical challenges [15] of the full NLO four-fermion cross section calculation. However, 
it should be mentioned that our calculation is restricted to the inclusive cross section, 
while a more flexible treatment of the final-state phase space is obviously desirable. This 
requires either applying effective-theory methods to four-fermion production amplitudes 
rather than the forward-scattering amplitude, or the consideration of specific cuts such as 
corresponding to invariant-mass distributions that allow for a semi-inclusive treatment. 
Interesting developments in this direction have recently been reported for top-quark pair 
production [48]. 
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A Coefficients of non-resonant contributions 

In this appendix we list the explicit expressions of the remaining coefficients C( h and K( 
in (37). The functions C( h are known analytically, and contain all the s-dependence of 
the cross section (except for the overall factor 1/s). They are determined by the photon 
and Z propagators and electroweak couplings. In the limit of vanishing fermion masses 
the only helicity configurations contributing to the cross section are h = LR,RL. The 
coefficients K{ are s-independent, and result from dimensionally regularized cut loop 
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integrals. Typically the last integration is performed numerically, after the subtraction 
of the singular terms which are integrated analytically, though some analytic results can 
be obtained. The results given below contain the contribution of the diagrams h4-h7 
in Figure 5 including their complex conjugates, except for cut h6, where the complex 
conjugate is the diagram itself, and cut h7, where the symmetric diagram is automatically 
taken into account by summing over the four flavours. 

Only the configuration e^e^ contributes to the cut diagram h4: 

r f - Q M 2 2 ( Qf , C e C f 

KJU = Kl = -0.266477 , K d M = K£ 4 = 0.190394 , (93) 



where Qf and Cf = WJ s " w - are the couplings of left-handed fermions to 7 and Z . Qf 
always denotes the charge of the particle (not the anti-particle) in units of the positron 
charge. For the cut diagram h5 we have 

CL„ - M$,st\-%+ C "°L + Cw Qf( ' ! ' '■ 



h:u > •"•"■"I S 2 ^ S ( s -M z )' s w s(s-M z ) s w {s-M z f 
K u h5 = K% = 0.455244 , K d h5 = = -0.455244 , (94) 

where C^ R = and Cf L = Cf = ~^Q e - In this case both left-handed and right- 
handed incoming fermions contribute (h = LR,RL), but only left-handed internal 
fermions. The coefficients of I16 are 



2 



K u m = < = Kl = K% = 0.0804075 , (95) 



while for h7 we get 



rf 0M*s(9f9l Q f c "°f QfCtcf c?c}c} 

^ h 7, h -™l w S w y g2 s{s _ M 2 z) s{s _ M 2) + (s _ M 2 )2 

Kl 7 = K d h7 = = K% = 0.0213082 , (96) 

where Qf, Qf and Cf, Cf are the couplings to 7 and Z of the particles in the same 
SU(2) doublet (i.e. fx, and u, d). 



B Hard one-loop coefficients 

We give here the explicit analytic results for the hard one-loop coefficients appearing in 
Section 4.1. 
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B.l Production vertices 

The general e~e + — > W~W + production operator we are concerned with in this appendix 
reads 

0, = ^C p (e 7 Me) (n*nj?) , (97) 



p M 2 



where C p = C Pi h is the hard matching coefficient and h = LR, RL refers to the helicity 
of the incoming leptons (e^e^ or e^ej). Starting with e^e R — > W~W + , the matching 
coefficient at tree level is equal to 1, as can be read off (13). At NLO we have 

C p ,lr = 1 + C£1 R + O (a 2 ) = 1 + ^c p % + O (a 2 ) , (98) 

where C ( p l \ R is the coefficient in (23). Before renormalization the NLO short-distance 
coefficient reads 

„(i),bare_ 1 ( 4M^y e 8cj + 10^ + 1 / AM 2 W 



(24 - 1) (244 + 164 - 1) M 2 W Co (0, M 2 ,, 0, Mf, M 2 W ) 



+ 



8cisf 



■'ww 

(24 - 1) Co (0, AM 2 ,, 0, 0, Mf, Mf ) 



2cis 



2 



_ ((4 + 174 - 16) M 2 H + Mjj M^ C (— Mfy, 0, 0, 0, Mgj 

4M£ 4 

(M| + Mgj M^ C (-M^, M^, 0, 0, M 2 H , M 2 ,) 
AM 2 H 4 

_ (24 + 324 + 324 - 114 - 16) M 2 W C (-M 2 ,, M 2 ,, 0, 0, Mf , M 2 W ) 

^ " 844 

3 (33 - 464) M 2 , Cp (M 2 ,, -M 2 ,, 0, 0, 0, Mjj 

84 

(44 - 1) (144 + 154 - 24 - 1) ^ Cp (Mjy, -M^, 0, 0, 0, M|) 

1644 

(1 - 24) 2 (4 + 1) (44 + l) 2 M^ Co (4M^, 0, 0, 0, 0, Mf) 



+ 



16<£s 



2 
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244 (1 - 44) s 



45 
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2444 

964-(10-24tt 2 )4-94-6 



24r 4 s 2 

it; uj 

(1284-644 + 44 + 234 + 5) m 



48c 4 ,s 2 
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where all functions appearing in the above expression, Co(p 2 ,p2,p 2 ,m 2 ,m?,,m 2 ) and 
£(g 2 , M 2 , M| ), are known analytically and are supplied in Appendix B.3. The coun- 
terterms in the scheme are computed from (46) and are given by 

c(1)iCt _ 44 - 224 " 1 ( 4M 2 , V e (Mjj - 3M 2 W M 2 H + 6M 4 ,) I (M 2 ,, M 2 H , M 2 W ) 



844 e V i' 1 ) i2.u lV 4 

(M| - bM 2 w ) t (0, Mg, M^) (84 + 274 - 5) I (0, Mf) 
12M 2 si 12c 2 - 2 



+ 



(424-114-1) £(M^,M^M|) 



1244 



_ (Mjj - AMyyMfj + 12M 4 ,) OBp (M 2 ,, M 2 ,, M 2 h ) 

24M^4 

+ (484 + 684 - 16 4 " 1) M 2 W dB {M 2 W) M 2 W , M 2 ) 



24c 4 , s 



2 



( 2M h ~ 3M 2 H M 2 W + 2M 4 ,) In M 4 3M 2 



— I i VLL. + 



24M*,(Af2 - M 2 w )s 2 w 12M 4 , s 2 w 16M 2 , s 2 w 

3m 2 (m 4 - Myy) In (l - ffi) 3^ 3m 4 
4M^ 4 8M 2 , 4 4M^4 

(124 - 724 + 264 -154 -2) In (jf) (44-224-1) In 2 
24r 4 s 4 4r 2 s 2 
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2(35 - 6in)cj + (-112 + 66t7r)4 + (13 + 3m)cj + 2 



W W 



The full renormalized coefficient is obtained by adding bare result and counterterms 

J 1 ) _„(!)> bare (l),ct /^-.n 

The poles of c^ LR are given explicitly in (51) and cancel once one takes into account soft 
and initial-state collinear radiation. 

Turning to the e R ej^ — > W~W + case, the matching coefficient C PjRL vanishes at tree 
level, as can be seen from (13). The NLO correction is therefore finite. We have 

C p ,rl = C { X + O (a 2 ) = ^ RL + O (a 2 ) , (102) 

where C^ RL is the coefficient in (23). We find 

(1) 44 Ml C (0, Mjy, -M^, 0, Mj, M 2 ,) _ g4 C (0,4M^,0,0, M|, M|) 
"■ WL 4(2^-1) 4(24-1) 

(244 + 204 - 5) 4 € (Mjy, M^, Mf) 2 (84 + 4 + 3) 4 i (AM 2 ,, Mf , Mf) 



c, 



+ 



34 (24 - 1) (44 - 1) 34(24-1) 



(644 + 44 + 1) 4 In ( ^ - l) (644 - 484 - 244 + 5) 4 In ^ 



v 1 ^ / + 1 w ~ w ~ w ' w v*v 



124 34 (24-1) (4c 

164 In 2 _ (324 + 44 + 1) 4 m 
3 124 

B.2 Virtual corrections to W decay 



(103) 



The decay of a W boson is implemented in the effective theory analogous to the produc- 
tion [9]. There are decay operators with collinear fields describing the decay products of 
the non-relativistic vector boson. For the flavour-specific decays under consideration we 
have up to NLO 

O d = -^== {Cm St /x C3 ,l7V C4 , l + C d:h n\u C3tLl % 4tL ) . (104) 

These operators would be needed for the calculation of the e _ e + — > /i~ v^ud scattering 
amplitude in the effective theory. However, for the total cross section (or the forward 
scattering amplitude) the directions c 3 , c 4 of the decay products will be integrated over 
and, as indicated in (4), there is no need to introduce collinear fields JJt C3 ,L, v c±,l-, u C3 ,l 
and d C4t L in the effective theory. The matching coefficients of the decay operators enter 
only indirectly through A^. The virtual correction to the W decay width is related to 
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the coefficient functions of the decay operators. Ignoring QCD corrections, at NLO we 
have 

C d , = 1 + C$ + O (a 2 ) = 1 + + O (a 2 ) , 

C d , h = 1 + C$ + (a 2 ) ee 1 + + O (a 2 ) . (105) 

We give here the explicit results for the electroweak corrections. The unrenormalized 
one-loop correction to the leptonic decay vertex reads 

Jl),ba*e 1 ( M^y 8<£ + 24 + 1 ( M 2 , 



4 + I) 2 (2c 2 w - 1) C (M 2 ,, 0, 0, 0, 0, Mf ) 



4 r 6 s 2 



+ (4 + 2) M^C p (M 2 ,, 0, 0, M 2 ,, Mf, 0) 



s 2 



+ 



(24 + 1) £ (Mw, Myyj Ml) « ~ 2c t + 1) ^ (jf ) 
2 s 2 " 4r 4 <? 2 

- (24 + 7T 2 ) c Q w + (tt 2 - 18 in) c% - 3mc 2 w + 6ivr + 6 



24cis 



2 



and the corresponding counterterms computed from (46) are 



-W'? / ^ x - e In ' ^ 



ji),ct_ c ^" 2^ + 1 / M'w\ c , 111 v mT; . 2c 2 +1 



Similarly the NLO bare correction to the hadronic vertex is given by 

(l),bare _J_ KV £ ,1K 

2e 2 V V? J 9e 2 V A* 2 



(l + 2c 2 )(l + 324) /M^V 6 , 1 ( M w X " 



+ 



72 44 e V ^ 2 / 3e V /i 2 

(84 + 184 + 114 - 1) Ml C (Mjy, 0, 0, 0, 0, Mf ) 

3644 

(4 + 2) M 2 , C (Ml, 0, 0, M 2 ,, Mf , 0) 



+ 

4 



+ (24 + 1 ) £ (M 2 ,, M 2 ,, Mf) (204 + 64 + 1) ^ 



(106) 



2 " 16c 2 s 2 e V ^ 2 ; 16c 2 4 3244 " 1 ; 



24 3644 
+ 1204 + (48 - 134tt 2 ) 4 - 6 + (244 + 224 + c* - 2) m ^ 
21644 ^ c w s w 
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and the corresponding counterterms are 



m,ct _ 164 - 50c* + 7 (M 2 W 



d ' h ' 2 1444 s 2 w e V V 2 

(164 - 324 + 7) In (|J) 16c 4 _ 50 4 + 7 



1444 4 28844 



(109) 



B.3 Integrals and auxiliary functions 



The results for the short- distance coefficients and their counterterms have been written 
such that all poles in e are apparent and the remaining functions are finite. We give 
here their analytic expressions. As usual the scalar two- and three-point functions are 
defined by 

B (k 2 m 2 m 2 ) " f [dl] M = ^ E ^ ddl h 10) 

B (k , mi ,m 2 ) _ J {p _ m 2 m + k y_ m 2y \ dl \ = ~^ > (HO) 

and 

[dl] 



C (kl,kl, (h + k 2 ) 2 ,m 2 1 ,ml,m 2 3 ) = J 



(Z 2 - m\){{l + - m 2 )((Z + k ± + k 2 ) 2 - mf) ' 

(HI) 

dBo(k 2 ,m 2 ,ml) is then defined as 

dB ^, mlml) ^_BM0^L\ q ^. (112) 

The auxiliary function £(k 2 , m 2 , m\) used in the expressions for the matching coefficients 
is related to the two-point function by 

B (k 2 ,mlm 2 2 ) = ±(^j + 2-t(k 2 ,m 2 1 ,m 2 2 ) (113) 

and satisfies £(/c 2 ,m 2 ,m 2 ) = £(/c 2 , m 2 ,, m 2 ) + \n(m\/m\). It is sufficient to give this 
function for the following special arguments: 



l(M 2 z ,M 2 w ,M 2 w ) = Ml 2M f W In (l + 



M 2 ZW - M\ 
2M 2 , 



+ Mz + Mzw ln ( 1 M 2 ZW + M 2 



2M 2 V 2M 2 , 

r , , r2 f2 , 2M 2 , -M 2 + M 2 ZW (M 2 - M 2 ZW 



m^M 2 w) ^-^ M2 ^ ^ 



2M| 



2M 2 W - M 2 - M 2 ZW ( M 2 + M 2 ZW 

2M 2 , V 2M 2 I ' 
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where we introduced M zw = a/M| — \.\ I 2 , \ /, 2 r . The explicit result for the derivative of 
the two-point function that is needed reads 

dB (M^,M^,M 2 z ) = (115) 

1 f M^-Mf ln /Mf\ Mf(3M^-Ml) ln / M 2 - Mf w 
M^l 2M^ \M^J M^M 2 ZW V 2M^M Z 

The analytic expressions of the finite three-point functions appearing in the results given 
in (99)-(109) can all be obtained from 

C a {0,Ml„-Ml;0,Ml,Ml,) = (116) 
1 ( 2 Li2 U - ?m + 2 Li 2 ( ffi-ffi - Li 2 (>£- 



M 2 / \ M 2 



C (0, 4M^, 0,0, M 2 Z , Mf) = (117) 

—L-i m 2 + m 2 f^M +2Ll2 fif 

SM 2 ^ V J V Ml J 2 \M ZW 

,AM 2 -M 2 \ fAM 2 -M 2 x 
+ 2Li 2 — ^ £ +2Li 2 — ^ ^ +7T 2 



/If 2 / V \/- 

JKZ W-.Z' / V 1V1 W+Z 



C (-M^, M^, 0, 0, M z , M^) = (118) 



1 1 lJ-^L^Vl,/ m - 



M 2 , - M 2 - M 2 ZW ) M^-Ml-M 2 ^ 
+ Li 2 < 5 ]-uJ M - ^ 



M 2 W - M 2 + M 2 ZW ) Z V M^-Mf + Mf^ 
C (M^,0,0,M^,M|,0) = (119) 



^\ V^f + ^iW V 2M 2 

2 !\/f2 n n n y\/f2^ 



VT 2 



C (M^,-M^,0,0,0,M 2 ) = (120) 
1 f , [2M 2 V , \ / /2M 2 , , ^ n . \ T . ( Mj 



In ^ + 1 In + 1 - 2m - Li 2 



4M 2 / [ V^I J\ \M 2 J J \ V ( 2M ^ + M I) 2 , 
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C (M^, 0,0, 0,0, Mf) = (121) 



where we introduced M w ± z = M w ± y/Myy — Mf . 
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